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ABSTRACT 


High  range  resolution  can  be  obtained  using  a  nuuibtr  of  methods  including  ultra- 
wideband  radar,  intrapulse  puise  compression,  and  super-resolution  techniques  This  thesis 
investigates  the  achievement  of  high  radar  range  resolution  by  means  of  the  step  frequency 
waveform  The  key  advantage  of  the  step  frequency  appioach  compared  to  other  methods 
is  the  wide  effective  bandwidth  resulting  in  high  range  resolution  while  still  maintaining  a 
narrow  instantaneous  bandwidth  which  eases  A/D  sampling  requirements  The  effects  of 
various  waveform  parameters  on  the  liigh  resolution  range  profile  are  detennined  for  both 
noise  and  clutter  limited  environments  A  novel  technique  is  developed  which  enhances  the 
resolution  of  moving  targets  in  clutter  and  estimates  the  target  velocity  This  method 
consists  of  clutter  cancellation  followed  by  successive  velocity  compensations  Finally, 
design  methods  are  developed  to  optimize  system  peiformance  in  a  low  PRF  mode  for 
both  noise  and  clutter  limited  environments 
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I.  INT».  ..  CTION 


One  important  feature  of  any  radar  system  is  the  resolution  of  the  target  in  the  range 
domain  A  fine  range  resolution  allows  the  user  to  distinguish  between  two  closely  spaced 
targets  that  might  otherwise  be  incorrectly  judged  as  a  single  target  and  can  aid  in  the 
target  identification  process  Moreover,  a  fine  range  resolution  will  enhance  the  target 
signal  relative  to  the  background  clutter  For  any  pi.ised  radar  system,  an  increase  in  the 
effective  waveform  bandwidth  will  improve  the  range  resolution 

In  recent  years,  a  number  of  high  resolution  techniques  have  been  available  to  the 
radar  designer  to  improve  the  range  resolution  of  the  radar  system  These  techniques  can 
be  grouped  into  four  main  categories:  ultra-wideband  radar  (UW'B),  superresolution 
techniques,  intrapulse  pulse  compression  techniques,  and  interpulse  pulse  compression 
techniques 

Ultra-Wideband  Radar  achieves  fine  range  resolution  because  of  the  exttemely 
narrow  pulses  (on  the  order  of  nanoseconds)  which  are  transmitted  without  carrier  As  a 
result,  the  large  instantaneous  bandwidth  imposes  severe  constraints  on  the  analog  to 
digital  conversion  process  since  the  waveform  sampling  must  meet  the  Nyquis;  criterion 
Moreover,  the  radar  requires  new  wideband  components  different  from  those  used 
in  conventio'.al  narrowband  radars  [Ref  I  p  31] 

Superresolution  techniques  involve  the  use  of  high-resolution  optimal  e.stimaiion 
methods  These  methods  are  based  on  an  eigenstruciure  technique  that  makes  use  of  the 
co\  ariance  matrix  of  the  received  signal  [Ref  2  p  18] 

Intiapulse  pulse  compression  (PC)  methods  have  been  commonly  used  to  increase  the 
range  resolution  of  radar  systems  These  techniques  increase  the  instantaneous  bandwidth 


by  applying  either  frequency  coding  or  phase  coding  to  each  transmitted  pulse  The 
received  pulses  are  processed  in  a  filter  matched  to  the  transmitted  signal  resulting  in  a 
compressed  pulse  Frequency  coding  techniques  include  the  common  linear  frequency 
modulation  (LFM)  or  chirp  method,  and  discrete  coded  segments  within  the  pulse  [Ref 
3:pp  420-424]  Phase  coding  techniques  include  binary  phase  codes  as  well  as  various 
polyphase  codes  [Ref  4:pp  223-228] 

The  step  frequency  waveform  is  an  interpulse  pulse  compression  method  in  which  the 
carrier  frequency  of  successive  pulses  is  varied  by  a  fixed  frequency  step  Wide  bandwidth 
is  indirectly  achieved  over  several  pulses  instead  of  within  a  single  pulse  as  in  conventional 
pulse  compiessicn  signals  The  effective  waveform  bandwidth  is  the  product  of  the 
number  of  coherently  integrated  pulses  N  and  the  frequency  step  size  AF 

The  reiurns  of  these  N  pulses  at  different  carrier  frequencies  constitute  the  frequency 
spectrutn  of  the  target  retlectiviiy  and  are  coherently  integrated  by  the  inverse  FFT  to  yield 
the  high  resolution  range  profile  (HRR  profile)  The  signal  processing  transforms  wide 
original  range  bins  into  finer  range  bins  The  KRR  profile  does  not  give  absolute  range 
position  of  the  target  but  gives  the  relative  position  of  each  scattering  center  within  the 
original  range  bin 

The  key  advantage  of  the  step  frequency  method  compared  to  other  pulse 
compression  techniques  is  that  the  range  resolution  is  increased  while  still  maintaining  a 
narrow  instantaneous  bandwidth  As  a  result,  the  analog  to  digital  conversion 
requirements  are  less  stnngent  Moreover,  as  opposed  to  the  UWB  waveform,  the  range 
resolution  can  be  achieved  by  merely  adding  a  stepped  frequency  synthesizer  to  a 
conventional  narrowband  radar  receiver  and  performing  th**  required  signal  processing  on 
the  received  signal 


The  step  frequency  w  iveform  has  two  primary  limitations  First,  as  opposed  to  a 
pulse-doppler  radar,  the  step  frequency  radar  cannot  directly  measure  the  velocity  of 
moving  targets  based  on  the  Doppler  shift  Second,  as  compared  to  intrapulse  pulse 
compression,  the  step  frequency  wavefonn  requires  additional  time  to  transmit  and  receive 
N  pulses  needed  to  obtain  the  equivalent  wide  bandwidth  of  a  single  narrow  pulse 

The  pur}2ose  of  this  thesis  is  to  investigate  the  means  by  which  high  range  resolution 
Can  be  achieved  via  the  step  frequency  waveform  Chapter  II  presents  the  primary  step 
frequency  radar  parameters,  a  brief  system  description,  and  the  signal  processing  that 
results  in  the  fine  range  resolution  Chapter  111  summarizes  the  relationships  between  the 
various  radar  design  parameters  and  tlieir  effects  on  the  overall  system  perlonnance 
Moreover,  this  chapter  presents  methods  for  designing  a  step  frequency  radar  system 
based  on  various  user  requirements  and  constraints  Chapter  IV  presents  the  results  of  a 
high  range  resolution  (HRR)  simulatioii  using  the  steo  frequency  waveform  The 
simulaiiOii  output  i‘,  iiir  puK.esscu  range  pronie  corresponding  to  one  original  range  bin 
The  HRR  simulation  graphically  displays  the  effect  of  varying  individual  radar  parameters 
in  both  noise-limited  and  clutter-limited  situations  In  addition,  a  new  technique  is 
presented  which  enhances  the  range  resolution  of  moving  targets  in  clutter  and  estimates 
the  target  velocity  by  clutter  cancellation  followed  by  repeated  velociry  compensations 
Appendix  A  contains  the  detailed  description  of  the  simulation  programs,  and  Appendix  B 
includes  the  program  codes 
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II.  OPERATING  PRINCIPLES  OF  STEP  FREQUENCY  RADAR 


A.  STEP  FREQUENCY  WAVEFORM 

The  step  frequency  waveform  can  be  described  as  an  inierpulse  version  of  the 
common  linear  FM  pulse  compression  waveform  A  series  of  N  coherent  pulses  are 
transmitted  whose  frequencies  are  monotonically  increased  from  pulse  to  pulse  by  a  fixed 
frequency  increment  AF  as  shown  in  Figure  2  1  Each  pulse  has  a  fixed  pulse  width  t,  and 
the  pulses  are  transmitted  at  a  fixed  pulse  repetition  frequency  (PRF)  [R'^f  5.pp  160- 
1 6 1 1  The  frequency  of  the  kth  pulse  is  given  by 

f,  =f,+(k-l)AF  (2  1) 


where  t"  is  the  nominal  carrier  frequency,  and  AF  is  the  frequency  step  size 

The  instantaneous  bandwidth  of  this  waveform  is  approximately  equal  to  the  inverse 
of  the  pulse  width  and  is  much  less  than  the  effective  bandwidth  The  wavefoim's  effective 
bandwidth,  denoted  as  in  contrast  to  the  instantaneous  bandwidth  is  determined  by 
the  product  of  the  number  of  pulses  N  and  ;  te  frequency  step  size  AF 

B,,  -  NAI  (2.2) 

The  range  resolution  for  any  waveferm  is  dependent  on  this  effective  bandwidth  of  the 
waveform  The  range  resolution  of  the  step  frequency  waveform,  in  meters,  is  gi\  en  as 
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Figure  2  I  Step  frequency  wavefonn  paranjeters  [Ref  6  o  23!' j 
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Ar  =  K  —  =K 

UbJ 


(23) 


where  K  is  a  constant  determined  by  the  weighting  window  and  is  greater  than  or  equal  to 
unity,  and  c  is  the  speed  of  light  (3  x  lO’'  m/s)  From  equation  (2  3),  the  range  resolution 
can  be  made  finer  by  either  increasing  the  number  of  pulses  and/or  increasing  the 
frequency  step  size  [Ref  6:p  234] 

The  coherent  processing  interval  (CPI)  is  the  total  time  duration  in  which  echo 
returns  are  collected  from  a  particular  target  for  signal  processing  The  CPI  is  simply  the 
product  of  the  number  of  pulses  N  and  the  pulse  repetition  interval  (PRI) 

CPI  =  N(PR1)  (2  4) 

The  coheient  processing  interval  is  also  referred  to  as  the  ramp  repetition  interval  (RRI)  or 
the  burst  time 

B.  SYS!  EM  DESCRIPTION 

The  step  frequency  radar  system  consists  of  elements  commonly  found  in  a  coherent 
pulse-doppler  radar  A  block  diagram  of  a  step  frequency  radar  system  is  shown  in  Figure 
(2.2)  A  pulse  is  generated  by  first  mixing  the  outputs  of  a  fixed  coherent  oscillator  at  the 
IF  frequency  range  (3U-60  MHz)  with  a  stepped  frequency  synthesizer  whose  frequency  is 
stepped  at  a  rate  equal  to  the  pulse  repetition  frequency  (PRF)  The  sum  frequency  from 
the  first  mixer  is  then  mixed  with  a  stable,  coherent  RF  frequency  oscillator  The  sum 
frequency  of  this  second  mixer  is  then  pulse  modulated  and  amplTied  prior  to  being 
transmitted 
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Figure  2  2  Block  diagram  of  step  frequency  radar  system  used  to  obtain  high 
range  resolution  (HRR)  target  profiles  [Ref  6  p  237] 
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The  carrier  frequency  of  each  transmitted  pulse  is  comprised  of  three  elements;  the  fixed 
IF  frequency  of  the  coherent  oscillator  (COHO),  the  fixed  RF  frequency  of  the  stable 
oscillator  (STALO),  and  the  variable  frequency  of  the  stepped  frequency  synthesizer 
Thus,  the  carrier  frequency  of  the  kth  pulse  is  given  by  the  following  equation 

fc,  (2  5) 

The  received  pulse  first  passes  through  the  duplexer  and  is  then  sent  through  an  RF 
amplifier  After  this  stage,  the  return  signal  is  mixed  with  the  STALO,  and  the 
downcom-erted  frequency  is  further  mixed  wiin  the  stepped  frequency  synthesizer  The 
stepped  frequency  synthesizer  is  synchronized  so  that  the  transmitter  and  receiver  are  on 
the  same  frequency  step  within  the  PRl  As  a  result  multiple  time  around  echoes  will  have 
frequencies  which  will  vary  by  multiples  of  the  frequency  step  size  AF  The  difference 
frequency  between  the  downconvened  STALO  output  and  the  second  luixer  (i  e  the 
stepped  synthesizer)  is  the  IF  frequency  The  signal  is  then  passed  through  an  IF  amplifier 
which  serves  as  a  matched  filter  for  the  echo  return  The  bandwidth  of  the  filter  is 
approximately  the  inverse  of  the  pulse  width  Thus,  multiple  time  around  echoes  w'hose  IF 
frequencies  fall  outside  the  filter  bandwidth  will  be  rejected  The  amplified  signal  is  then 
split  into  two  separate  channels  The  first  channel  is  mixed  with  the  coherent  oscillator 
directly,  and  the  second  channel  is  mixed  with  the  COHO  signal  after  it  has  undergone  a 
90  degree  phase  shift  The  difference  frequencies  in  channels  one  and  two  comprise  the  in- 
phase  (1)  and  quadrature  (Q)  components  respectively  Each  of  the  1  and  Q  components 
are  in  the  video  frequency  range  [Ref  6  pp  236-237] 

Both  the  in-phase  and  quadrature  signals  are  sent  through  an  analog  to  digital 
converter  (A.'D)  whose  outputs  are  sampled  at  a  rate  equal  to  the  inverse  of  one  pulse 


width  The  I  and  Q  samples  for  each  range  bin  are  stored  in  memory  until  all  the  pulses 
witliin  the  burst  h  ve  been  received  for  processing  Complex  samples  for  each  ran  ;  bin 
are  transformed  by  an  inverse  FFT  to  form  the  high  resolution  range  profile  (HRR  profile) 
for  that  range  bin  To  reduce  sidelobes  in  the  range  profile,  the  I  and  Q  samples  are  usually 
windowed  prior  to  computing  the  inverse  FFT 

One  of  the  key  advantages  of  a  step  frequency  radar  compared  to  radars  using  other 
wideband  waveforms  is  that  the  naiToxv  instantaneous  bandwidth  eliminates  the 
requirement  for  extremely  hi  'h  analog  to  digi.^l  conversion  rates  which  can  be  a  major 
limiting  factor  in  the  system  design  The  major  limitation  in  the  step  fiequency  method  is 
the  time  required  to  transmit  and  receive  N  pulses  to  achieve  the  desired  fine  range 
resolution 

C.  W'AVTFORM  PROCESSING 

As  opposed  to  conventional  dars,  the  pulses  within  one  coherent  processing  interval 
have  different  carrier  frequencies  The  target  return  from  a  pulse  at  a  specific  frequency  is 
the  magnitude  of  the  target  reflectivity  at  that  frequency  The  collection  of  these  samples 
from  N  pulses  at  N  frt  ^uencies  essentially  constitutes  a  frequency  rpectrum 

The  return  signal  consists  of  a  video  signal  comprised  of  in  phase  and  quadrature 
components  as  given  by 


I(k)  =  cosd>^ 

(26) 

Q(k)  =  A(^  sin 

(2  7) 

0 


and  the  complex  samples  are  then  formed  as 


S(k)  =  I(k)-jQ(k)  =  A,exp(-j(DJ 


(2.8) 


where 


=  27rf,t, 


(2  9) 


(2.10) 


and  f,  refers  to  the  kih  frequency  step  of  the  frequency  stepped  synthesizer,  and  t^ 
denotes  the  round  trip  time  due  to  a  stationary  scatterer  at  range  R.  Processing  of  received 
signals  can  be  interpreted  as  converting  frequency  domain  target  reflectivity  data  via  the 
inverse  Fourier  transform  into  the  time  domain  to  form  the  high  range  resolution  (HRR) 
profile 

If  the  target's  velocity  is  in  the  radar's  line  of  sight,  the  target’s  range  can  be 
generalized  to  become 


R,  =  R„-t  (k-l)vT  (2.11) 

where  refers  to  the  initial  range  at  which  the  target  is  detected,  v  is  the  radial  velocity 
of  the  target,  and  T  equals  the  PRI  A  positive  radial  velocity  indicates  an  opening  target. 
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and  a  negative  radiai  velocity  indicates  a  closing  target  Moreover,  a  single  point  scatterer 
with  a  constant  velocity  is  assumed  in  these  equations 

The  overall  expression  for  the  baseband  return  signal  can  be  determined  by 
substituting  the  instantaneous  frequency  from  equation  (2.!),  the  phase  from  equation 
(2  9),  the  round  trip  time  from  equation  (2.10),  and  the  target  range  from  equation  (2  11) 
into  the  expression  for  s(k)  in  equatioti  (2  8)  The  resulting  expression  is  given  by 

S(k)  =  A,  exp(-j  — [fXk  -  1)AF][R,  +(k  -  l)vT])  (2  12) 

c 

Expanding  the  expression  inside  the  exponent  in  equation  (2  12)  results  in  a  nonlinear 
phase  term  v/hich  causes  range-doppler  coupling  for  moving  targets  The  degree  of  range- 
doppler  coupling  profoundly  impacts  the  range  resolution  of  moving  targets,  as  well  as  the 
integration  gain  of  the  signal  processing  [Ref  7:pp  4-5]  A  complete  description  of  these 
effects  will  be  discussed  in  a  later  section  The  set  of  N  sampled  baseoand  signals,  S(k), 
from  a  given  range  bin  of  width  ct/2  are  converted  to  a  time  domain  range  profile  by  the 
inverse  FFT  [Ref  6:p  238] 


s(n)  =  ^^S(k)exp|^j^^^j,  0<n<N-l  (213) 

1.  Phase  Relationships  of  Return  Signals  of  Successive  Pulses 

The  phase  of  a  target  return  from  a  given  range  for  a  specific  pulse  is  obtained 
from  equations  (2  9)  and  (2  10)  and  is  given  by 
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0),= 


(2  14) 


47if,R, 
c 

The  step  fiequency  method  produces  an  induced  phase  shift  between  successive  pulses 
because  of  the  shift  in  the  carrier  frequency  from  pulse  to  pulse.  The  time  rate  of  change  of 
the  phase  is  given  by: 


d(<I),)  ^  47td(f,R,)  ^  4nR,  ^ 

dt  c  dt  0*^0“  (215) 


The  corresponding  frequency  is  given  as 


f .  = 


2rt  dt  c  ‘‘  c 


(2  16) 


The  rate  at  which  the  phase  shifts  from  pulse  to  pulse  depends  on  two 
components,  as  in  equation  (2  16)  The  first  component  is  the  rate  of  change  of  the 
instantaneous  carrier  frequency,  f^  The  second  component  is  the  Doppler  shift  of  the 
target  caused  by  relative  motion  in  the  radial  line  of  sight  If  the  relative  motion  between 
the  target  and  the  radar  is  negligible,  equation  (2  16)  can  be  approximated  as  follows 
[Ref  6:p  239] 


^  ^  2R,  AF  2R,  ^ 

f.  =  f. 

c  T  c 


PRf 


(2  17) 
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The  maximum  unambiguous  range  is  that  for  which  the  induced  phase  shift 
frequency  equals  the  PRF  Therefore,  the  maximum  unanibig»ious  range  of  the  high 
resolution  range  profile  is  given  a'^: 


R 


u 


c 

2^ 


(2  18) 


The  unambiguous  range  is  inversely  proportional  to  the  frequency  step  size  The 
above  expression  is  similar  to  the  unambiguous  range  for  a  conventional  radar  except  that 
the  frequency  step  size  appears  in  the  denominator  instead  of  the  PRF.  If  is  1  MHz,  the 
unambiguous  range  will  be  150  meters  which  is  very  small  in  comparison  to  conventional 
radars  As  a  consequence,  any  target  with  a  range  which  is  a  multiple  of  the  unambiguous 
range  will  be  folded  into  the  range  profile  resulting  in  range  amb’guity  For  example, 
targets  with  ranges  of  100  meters,  250  meters.  400  meters,  etc  will  a!!  appear  in  the  same 
location  in  the  HRR  profile  Therefore,  one  cannot  obtain  absolute  ranges  to  individual 
scatterers  Instead,  the  HRR  profile  gives  relative  range  positions  of  scatteiei::,  contained 
within  one  pulse  width  [Ref  6:p  243] 

2.  Range  Bin  Transformation 

The  original  range  resolution  cell  of  the  step  frequency  radar  prior  to  taking  l  : 
inverse  FFT  is  determined  by  the  pulse  width  and  is  given  by 


Ar, 


cx 

T 


(2  19) 


where  t  is  the  pulse  width  (if  pulse  compression  has  been  used  earlier,  x  refers  to  the 
compressed  pulse  width)  As  a  result  of  the  signal  processing,  the  unambiguous  range 
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window  of  width  is  split  into  N  fine  range  resolution  cells  whose  width  (assuming  no 
range  weighting  functions  so  that,  K,  the  tesolution  loss  factor  due  to  window. ng,  is  set  to 
unity)  is  equal  to 


N  2NAF 


(:.20) 


The  fraction  of  the  unambiguous  range  window  which  is  occupied  by  the  original 
range  bin  Ar^^  is  given  as  follows 


K 


cx 

■2 

c  ' 
2AF 


=  xAF 


(2  21) 


The  original  range  bin  is  subdivided  into  finer  range  bins  The  number  of  processed  range 
bins  corresponding  to  the  original  range  bin  is  equal  to 


-^  =  (tAF)-N 
Ar 


(2.22) 


3.  The  Implications  of  the  Product  tAF 

The  product  of  the  pulse  width  and  the  frequency  step  size  has  major 
consequences  on  the  nature  of  the  HRR  profile  In  general,  there  are  three  possible 
outcomes 

When  tAF  is  equal  to  one,  the  unambiguous  range  window,  R  ,  equals  the 
original  range  bin  size  of  cx/2  In  effect,  the  old  range  bin  is  split  into  N  segments  with 
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each  segment  having  a  range  resolution  equal  to  Ar  This  situation  represents  the  minimum 
required  sampling  rate  in  the  frequency  domain  for  a  given  pulse  width  It  should  be  noted 
that  the  sampling  rate  in  the  frequency  domain  is  equal  to  the  frequency  step  size  The 
drawback  to  setting  tAF  to  unity  is  that  it  becomes  difficult  to  notice  the  wraparound  of 
scatterer  peaks  caused  by  slight  changes  in  target  range  Moreover,  aliasing  may  occur  if 
there  is  interference  from  adjacent  original  range  bins  [Ref  7  ;p  7]  Finally,  if  the  original 
range  cell  contains  a  moving  target  and  clutter,  the  target  will  be  masked  by  clutter  [Ref 
5pp  168-169] 

When  tAF  is  less  than  one,  the  original  range  bin  comprises  only  a  fraction  of 
the  unambiguous  range  window  This  situation  represents  oversampling  in  the  frequency 
domain  and  leads  to  a  loss  of  resolution  in  the  range  profile.  However,  the  unused  portion 
of  the  range  profile  repicsents  a  clutter  free  zone  which  can  be  used  to  detect  moving 
targets 

When  tAF  is  greater  than  one,  the  unambiguous  range  window  is  smaller  than 
the  original  range  bin  This  situation  represents  undersampling  in  the  frequency  domain 
As  a  result,  those  regions  of  the  original  range  bin  not  contained  within  the  unambiguous 
range  w  indow  w'ill  fold  over  resulting  in  aliasing  of  the  HRR  profile  This  situation  must 
be  avoided  for  mapping  range  profiles  [Ref  8:pp  602-603] 

One  basic  requirement  in  choosing  the  freauency  step  size  is  that  the 
unambiguous  range  w'indow,  R^,  must  encompass  the  maximum  target  extent  E  to  avoid 
target  wraparound,  as  shown  below,  or  expressed  alternatively  in  equation  (2  24) 


(2  23) 
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(224) 


AF  < 


c 

2E 


An  additional  requirement  for  good  detectability  is  that  the  pulse  width  must 
encompass  the  maximum  target  extent  If  the  pulse  width  does  not  encompass  the  target 
extent,  different  scatterers  within  the  target  will  not  be  located  within  the  same  original 
range  bin  resulting  in  a  loss  of  SNR  [Ref  6  p  245]  To  encompass  the  target  extent,  the 
pulse  width  must  be  greater  than  or  equal  to  the  round  trip  time  from  the  nearest  point  of 
the  target  to  the  farthest  point  on  the  target,  given  by 

2F 

T^—  (2  25) 

c 

D.  EFFECTS  OF  RADIAL  VELOCITY  ON  THE  RANGE  PROFILE 

Radial  velocity  between  the  radar  and  the  target  results  in  additional  pulse  to  pulse 
phase  shift  caused  by  the  Doppler  effect  The  second  component  in  equation  (2  16) 
represents  the  Doppler  frequency  shift  and  can  be  rewritten  as 


(226) 


where  is  the  wavelength  of  the  kth  pulse,  and  v  is  the  relative  radial  velocity  between 
the  radar  and  the  target  The  additional  phase  shift  due  to  the  radial  velocity  results  in  a 
mismatch  in  the  inverse  FFT  which  acts  as  matched  filter  for  the  stepped  frequency 
waveform  only  if  the  radial  velocity  is  equal  to  zero  [Ref  7  p  44]  The  consequences  of 
the  mismatch  in  the  inverse  FFT  have  three  major  effects  on  the  range  profile. 
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The  first  effect  is  an  attenuation  of  the  peak  signal  level  of  the  scatterer  Since  the 
mean  noise  level  remains  the  same,  this  results  in  a  reduction  in  the  signal  to  noise  ratio 
(SNR)  The  second  effect  is  the  dispersion  or  spreading  of  the  peak  signal  into  adjacent 
range  bins  The  amount  of  dispersion  is  determined  by  a  dimensionless  parameter,  P, 
which  represents  the  number  of  range  bins  the  signal  spreads  in  one  coherent  processing 
interval: 


P^vN(TR^  (2  27) 

Ar 

The  effects  on  the  signal  peak  for  different  values  of  P  are  shown  in  Figure  2  3  In  general, 
the  dispersion  factor,  P,  should  be  kept  below  three  for  reasonable  performance  [Ref  7  p 
5]  In  effect,  the  dispersion  of  the  signal  results  in  a  loss  of  range  resolution  because  the 
peak  has  broadened  from  one  range  bin  of  width  Ar  to  several  range  bins  whose  width  is 
equal  to  PAi  One  should  note  that  the  dispersion  can  be  effectively  reduced  by  increasing 
the  PRF  of  the  radar,  without  having  to  adjust  the  number  of  pulses  N  or  the  frequency 
increment  AF  [Ref  7  p  46] 

The  third  effect  of  radial  velocity  is  a  shifting  of  signal  peaks  ii.  the  HRR  profile  This 
effect  is  similar  to  the  range-doppler  coupling  effects  in  conventional  linear  FM  pulse 
compression  waveforms  In  addition  to  the  dispersion  or  broadening  of  the  peak,  the 
location  of  the  peak  is  shifted  by  L  range  bins  where  L  represents  the  bin  shift  factor  and  is 
given  by  [Ref  7  p  45] 
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(228) 


L  = 


f.N(PRI) 


where 


r  =f.  +4B 


«n 


(229) 


Thi  number  of  shifted  rang£  bins,  L,  is  directly  proportional  to  the  dispersion  factor 
P  The  factor  is  the  median  va>’ic  of  the  carrier  ftequency  within  one  coherent 
processing  interv'al  [Ref  7;p  45],  If  the  sum  of  the  bin  shift  factor  and  thi.  original  bin 
location  of  the  signal  peak  exceeds  the  N  pulses  usrd  in  the  inverse  FFT,  the  signal  peak 
will  wraparound  from  the  beginning  of  the  HRR  resulting  in  an  aliased  piofile.  It  will  be 
difficult  to  te'l  from  the  riRR  whether  the  target  peak  is  in  its  actual  position,  or  has  been 
wrapped  around  several  times  due  to  radial  motion  The  bin  shift  factor  represents  the 
number  of  shifted  range  bins  from  the  original  peak  location  to  the  beginning  of  the 
broadened  shifted  peak,  not  the  ce..ter  of  the  shifted  peak 

One  Avay  to  reduce  the  effects  of  radial  velocity  is  to  apply  velocity  compensation  to 
the  frequency  domain  data  prior  to  taking  the  inverse  FFT.  Assuming  that  all  scatterers 
within  a  single  target  have  the  same  known  rad  a!  velocity,  one  can  cancel  the*  mismatch  by 
rr.ultiplying  the  returns  front  i.'ie  appioximaie  range  location  of  the  target  by  a  correction 
factor  iRef  6p  253]  The  expressiofi  lor  the  correction  factor  is  given  as 


V  (k)  =  exp(i^v(f  ^{k-l),4FHk-l)T)  (2  30) 

c 
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where  V^(k)  is  the  compensation  applied  to  the  signal  return  of  the  kth  pulse  by 
multiplying  it  with  S(k)  The  key  limitation  in  applying  a  compensation  factor  is  that  the 
target  velocity  must  be  knowm  in  advance  by  another  method,  or  be  given  by  another 
radar  Errors  in  the  compensated  speed  are  most  significant  if  the  sign  of  the  velocity  is 
estimated  incorrectly  For  example,  if  the  true  target  velocity  is  +300  m/sec  (opening 
target)  and  the  compensated  velocity  is  -300  m/sec  (closing  taiget),  the  resulting 
dispersion  will  prevent  the  target  peak  from  being  observed  in  the  range  profile 

In  summary,  the  step  frequency  waveform  consists  of  transmitting  and  receiving  a 
series  of  N  coherent  pulses  whose  frequencies  change  by  a  fixed  amount  AF  The  wide 
effective  bandw'idth  equals  the  product  of  the  number  of  pulses,  N,  and  the  frequency  step 
size,  AF  It  is  this  effective  bandwidth,  achieved  over  several  pulses,  which  results  in  the 
fine  processed  range  resolution  The  step  frequency  radar  system  is  similar  to  conventional 
narrowband  systems  with  the  addition  of  a  stepped  frequency  synthesizer 

Frequency  domain  reflectivity  data  are  coherently  processed  by  the  inverse  FFT  to 
produce  the  high  resolution  range  profile  for  each  original  range  bin  The  fraction  jf  the 
HRR  profile  occupied  by  the  original  range  bin  is  gis'en  by  the  product  tAF  This  product 
m.usi  be  less  than  or  equal  to  unity  to  prevent  aliasing 

Radial  velocity  between  the  radar  and  the  target  results  in  dispersion  of  the  target 
range  profile,  target  range  shut  caused  by  range-doppler  coupling,  and  a  reduction  in 
SNR  Proper  velocity  compensation  corrects  these  deficiencies  but  requires  advance 
knowledge  of  the  correct  target  velocity  If  velocity  compensation  is  not  possible, 
increasing  the  PRF  will  alleviate  the  effects  of  target  motion  on  the  IIRR  profile 
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m.  DESIGN  METHODOLOGY  FOR  STEP  FREQUENCY  RADAR 


The  design  of  a  step  frequency  radar  system  requires  a  thorough  understanding  of  the 
efiects  of  the  primary  design  parameters  on  the  overall  system  performance  These 
primary  design  parameters  consist  of  the  following 

1  N  ,  the  number  of  pulses  witliin  a  coherent  processing  interval 

2  AF,  the  frequency  step  size 

3  fp^,,  the  pulse  repetition  frequency 

4  T,  the  pulse  width  (compressed  pulse  width  if  pulse  compression  is  used) 

Once  »he  effects  of  these  individual  parameters  are  known,  then  it  is  possible  to  determine 
optimal  combinations  of  parameters  to  achieve  certain  specified  design  goals  The  design 
process  is  shaped  by  the  operational  requirements  determined  by  the  user,  as  well  as  the 
environment  in  which  the  radar  will  function 

A.  EFFECTS  OF  INDIVIDUAL  PARAMETERS 

In  this  section,  a'l  other  parameters  are  assumed  to  be  held  constant  while  the 
parameter  of  interest  is  being  increased 

1.  Increase  N,  the  number  cf  pulses 

An  increase  in  the  number  of  frequency  steps  N  within  a  coherent  processing 
interval  will  increase  the  effective  wavefonn  bandwidth 

B,,  =  NAF  (3  1 ) 


21 


The  increased  bandwidth  will  result  in  a  finer  processed  range  resolution,  Ar,  since  the 
range  resolution  is  inversely  proportional  to  the  waveform  bandwidth: 


At  = 


Kc 


(3.2) 


An  increase  in  N  will  also  lead  to  a  longer  col'.erent  processing  interval  (CPI),  which  is  the 
product  of  the  number  of  pulses  in  a  burst  and  the  PRJ  which  requires  a  longer  dwell  time 
The  longer  processing  time  will  also  lead  to  a  finer  frequency  resolution  in  the  FFT  since 

Af  (3  3) 

N 

I  or  moving  targets,  an  increase  in  the  number  of  pulses  will  lead  to  a  lar  ge  increase  in  the 
dispersion  factor  P  (  P  x  )  and  an  increased  bin  shift  factor  L,  given  by 


P  =  = 


vN 


2NAF 


2vN-AF 

fp*.pKc 


(34) 


L  = 


vN 


f 


vN 


2f  vN 


NAF.  _K^  f,^Kc 


2NAF 


(3  5) 
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2.  Increase  AF,  the  frequency  step  size 

At.  increase  in  the  frequency  increment  between  pulses  will  increase  the 
waveform  bandwidth  as  in  equation  (3  1)  and  results  in  a  finer  processed  range  resolution 
as  in  equation  (3.2)  In  addition,  increasing  the  frequency  step  size  will  decrease  the 
unambiguous  range  window  since 


R  = 


2AF 


(3  6) 


For  nsoving  targets,  increasing  the  frequency  step  size  will  increase  the  dupersion  of  the 
target  peak  as  in  equation  (3  4)  However,  AF  does  not  influence  the  bin  shift  factor  L  as 
indicated  in  equation  (3  5)  Increasing  the  frequency  spacing  between  pulses  without 
regard  for  the  pulse  width  can  lead  to  aliasing  in  the  range  profile  if  the  product  of  tAF 
exceeds  unity.  Aliasing  will  also  occur  if 


(3  7) 


3.  Increase  the  PRF 

An  increase  in  the  PRF  will  lead  to  both  improved  and  earlier  detection  of 
targets  The  radar  range  equation  written  in  terms  of  the  peak  power.  P,,  is  given  bv 


R^ 


P.G  X'na 
(47t)’kT,B,FL 


(3  8) 
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where 


P,  =  the  peak  power  of  the  transmitter  in  watts 

G  =  the  antenna  gain  of  the  transmit/receive  antenna  (assumed  to  be  equal) 
X  =  the  radar  wavelength  in  meters 

n  -  the  number  of  pulses  coherently  integrated  within  one  scan 
a  =  the  radar  cross  section  of  the  target  in  square  meters 

R  =  the  detection  range  of  the  target  in  meters 

k  =  Boltzmann's  constant  (1.38  xlO'*^) 

=  the  noise  temperature  in  degrees  Kelvin 
Bf  =  the  radar's  instantaneous  bandwidth  in  Hertz 
F  =  the  receiver  noise  figure 

L  =  a  loss  factor  incorporating  all  other  system  losses 

Equation  (3  8)  can  be  rewritten  as: 


V'nL„  P/„,„(47r)’kT„FLfp^ 


(3  9) 


where 


(3  10) 
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As  indicated  in  equation  (3.9),  the  minimum  SNR  needed  to  detect  a  target  at  a  given 
range  R,  is  inversely  proportional  to  the  PRF  resulting  in  improved  detection  as  the  PRF  is 
increased 

Increasing  the  PRF  will  also  reduce  the  range*doppler  coupling  for  moving 
targets  since  both  the  dispersion  and  bin  shift  factors  are  inversely  proportional  to  the  PRF 
as  shown  in  equations  (3.4)  and  (3  5)  For  noise  limited  situations,  the  reduction  in  range- 
doppler  coupling  eases  the  velocity  compensation  requirements  and  reduces  the  chances  of 
aliasing  for  uncompensated  range  profiles  Although  increasing  the  PRF  will  reduce  the 
chances  of  aliasing  in  clutter  limited  situations,  it  may,  however,  prevent  the  migration  of 
the  target  peaks  into  the  clutter  free  region  v/hich  is  required  for  target  detection  in  the 
presence  of  strong  clutter  The  PRF  has  no  effect  on  the  range  profiles  of  stationai'y 
targets  in  either  a  noise  limited  or  clutter  limited  situation  Furthermore,  the  PPJ  has  no 
effect  on  the  processed  range  resolution  Ar  or  the  unambiguous  range  window  as  in 
equations  (3.2)  and  (3.6)  respectively  The  inverse  FFT  operation  yields  the  HRR  profile 
containing  N  range  bins  with  a  resolution  of  Ar,  which  does  not  depend  on  the  PRF 
4.  Increase  t,  the  pulse  width 

An  increase  in  the  pulse  width  will  increase  the  energy  or  average  power 
directed  at  the  target  which  in  turn  will  lead  to  both  earlier  detection  for  a  fixed  SNR,  or  a 
reduced  minimum  SNR  to  detect  a  target  at  a  given  range,  as  follows 


R" 


(P,T)G'X’no 


EG'X'no 


(4n)^kT,FL  {4a)-'kT,FL 


(3.11) 
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(3  12) 


xz 


2E 


c 


However,  for  fine  range  profiling  (i  e  mapping  of  stationary  targets),  this  condition  need 
not  be  met 


B.  SUMMARY  OF  EQUATIONS  RELATING  DESIGN  PARAMETERS 

Equations  relating  the  important  design  parameters  are  summarized  below  in  Table 
3.1. 


TABLE  3  1  EQUATIONS  RELATING 


Design  Parameter 


Processed  Range 
Resolution 


Equation 


Effective  Waveform 
Bandwidth 


Unambiguous 
Range  Window 


Ar 


Kc 

2NAF 


B.,  =  NAF 


R,  = 


2AF 


MPORTANT  DESIGN  PARAMETERS 


Design  Parameter 


Dispersion 

Factor 


Range  Bin  Shift 
Factor 


Clutter  Free 
Fraction  of  Range 
Profile 


Equation 


P  = 


vN 

fpRf 


R 


=  xAF 


I 

i 

I 

I 

I 

I 

i 

j 
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C.  DESIGN  METHODS  FOR  THE  STEP  FREQUENCY  WAVEFORM 


The  following  section  contains  various  design  metiiods  to  achieve  high  range 
resolution  in  either  noise  limited  or  clutter  limited  situations  using  the  step  frequency 
waveform  For  each  method,  the  principal  requirements,  the  adjustable  parameters,  and 
key  assumptions  will  be  stated  prior  to  the  detailed  description  of  the  various  design  steps 
Each  of  these  methods  assumes  a  low  PRF  mode  of  operation,  which  results  in 
unambiguous  initial  target  ranges  out  to  the  maximum  detection  range. 

1.  Design  Method  No.  i:  clutter  environment 

This  design  method  is  based  on  the  assumption  that  clutter  is  the  limiting  factor 
in  the  detection  of  targets  The  effects  of  thermal  noise  on  the  range  profile  are  minimal 
compared  to  clutter  This  situation  would  be  expected  if  the  main  beam  is  pointed 
downward  The  principal  requirements  are  as  follows 

1  Processed  range  resolution,  Ar  (meters) 

2  Maximum  target  detection  range  (unambiguous),  (meters) 

Maximum  target  extent,  E  (meters) 

Minimum  and  maximum  target  radial  velocities,  v„„„  (m/s) 

The  adjustable  design  parameter.;  are  as  follows 

1  Frequency  step  size,  AF  (Hz) 

2  Pulse  width,  x  (seconds) 

3  Number  of  pulses  within  coherent  processing  interval,  N 

4  Pulse  repetition  frequency,  PRF  (Hz) 

5  Clutter  Free  portion  of  HRR  (expressed  as  fraction  of  entire  profile) 

6  Time  on  Target,  T,^.  (seconds) 
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The  design  process  consists  of  the  following  steps  in  the  order  shown  below 

1  From  the  required  processed  range  resolution,  calculate  the  waveform  bandwidth; 

B..  =  NaF  =  ^  (3  13) 

2  Determine  the  pulse  width  necessary  for  the  pulse  to  encompass  the  target  for  good 
detectability  However,  a  larger  pulse  width  will  also  increase  the  amount  of  clutter 
in  the  range  cell 


t 


c 


(3  14) 


3  Determine  the  percentage  of  clutter  free  space  in  the  range  profile  The  fraction 
of  the  range  profile  containing  clutter  is  specified  by  the  following  equation 


ct 

-^  =  (tAF)  (3  15) 

Choosing  a  low  value  for  the  product  xAF  (i  e  0  1)  will  increase  the  fraction  of 
clutter  free  space  which  will  result  in  improved  detection  of  moving  targets 
However,  it  will  worsen  the  processed  range  resolution  Conversely,  choosing 
values  of  xAF  close  to  unity  will  eliminate  the  clutter  free  space  where  targets 
may  be  detected  Also,  fast  moving  targets  are  likely  to  wrap  around  for  high  values 
of  tAF,  which  will  degrade  their  detectability  The  product  xAF  must  not  exceed 
unity  in  order  to  prevent  aliasing 
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4  The  frequency  step  size  can  be  specified  as  the  following: 


fC*id 


(3  16) 


The  value  of  AF  determined  by  equation  (^16)  must  also  satisfy  the  following 
requirement  shown  below  to  prevent  aliasing 


AF< 


c 


(3  17) 


If  the  value  of  the  step  size  exceeds  the  constraint  of  equation  (3  1 7),  the  appropriate 
course  of  action  is  to  reduce  the  frequency  step  size  so  that  equation  (3,1 7)  is 
satisfied  and  then  to  increase  the  pulse  width  as  needed  to  maintain  the  desired 
amount  of  clutter  liec  space  (i  e  keep  tAF  constant) 

5  Determine  the  number  of  frequency  steps  N  based  on  the  required  bandwidth 
and  the  computed  frequency  step  size 

(3.18) 

AF 


6.  Select  the  PRF  based  on  the  following  considerations 

A  The  upper  bound  will  be  determined  by  the  PRF  needed  to  move  the  slowest 
target  from  the  beginning  of  the  clutter  region  to  just  outside  the  clutter  region 
i  T  he  number  of  bins  in  clutter  region  equals  (tAF)N 
(one  pulse  width) 
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ii  The  required  amount  of  bins  to  be  shifted  equals  (tAF)N 
(to  move  the  target  outside  the  clutter  legion) 
iii.  The  bin  shift  factor  L  is  given  as: 


L  = 


fcVm.,N 


(3.19) 


XJf^ 


iv  Set  the  bin  shift  factor  equal  to  the  number  of  bins  to  be  shifted: 


f  V  N 

- ^  (tAF)N 


)uppei 


(3.20) 


V  After  simplification,  the  upper  bound  on  the  PRF  is 


ftv. 


—  c  mtft 


BjjfArAFt 


(3.21) 


B  The  lower  bound  on  the  PRF  will  be  determined  as  that  PRF  which  will  cause 
the  target  at  the  end  of  the  clutter  region  with  the  maximum  expected  velocity 
to  the  upper  edge  of  the  unambiguous  range  window  without  wrap  around 

i  The  number  of  range  bins  from  the  end  o^the  clutter  region  to  the 
upper  edge  of  the  range  profile  equals  N(  1  -  tAF) 

ii  Set  this  above  quantity  equal  to  the  expression  for  the  bin  shift 
factor  L,  and  solve  for  the  resulting  PRF 
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(3.22) 


f  V 

ff  _ 

V  PRf  f  nun  ^  A  y  1  A  r*  \ 

B^Ar(l-TAF) 

If  (tAF)  »  0.5,  and  target  velocity  is  constant,  the  upper  and 
lower  bounds  are  equal 
The  underlying  rationale  in  determining  the  upper  and  lower  bounds  of  the  PRF  is 
that  an  increase  in  the  PRF  will  reduce  the  range-doppler  coupling  resulting  in  a 
reduction  in  both  the  number  of  bins  that  the  target  peak  shifts  and  the  dispersion  of 
the  peak  As  a  result,  the  PPJ^  must  not  be  so  high  to  the  point  where  the  bin  shift  is 
insufficient  to  cause  the  moving  target  to  migrate  to  the  clutter  tree  zone 
Conversely,  too  low  a  PRF  will  increase  the  range-doppler  coupling  which  nay  cause 
the  target  peak  to  wrap  around 

7  The  time  on  target  is  determined  by  the  PRF  and  the  number  of  frequency  steps 


T 


(3  23) 


The  PRF  that  is  chosen  must  lie  in  the  range  specified  by  equations  (3  21)  and 
(3  22)  However,  the  value  chosen  must  result  in  an  unambiguous  range  greater  than  or 
equal  to  the  maximum  target  detection  range  If  no  PRF  value  within  the  range  satisfies 
the  unambiguous  maximum  target  range  requirement,  then  the  system  must  be  redesigned 
for  a  different  range  of  target  velocities  or  for  a  different  maximum  target  detection  range 
2.  Design  Method  No.  2:  noise  environii.ent 

This  design  method  is  based  on  the  assumption  that  the  effects  of  clutter  are 
negligible  This  assumption  is  valid  in  situations  where  the  main  beam  is  elevated  upward 
The  design  steps  do  not  assume  that  velocity  compensation  has  been  applied  to  the  profile 
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Velocity  compensation  can  be  included  at  the  discretion  of  the  radar  designer  if  the 
performance  requirements  cannot  be  me;  otherwise  The  principal  requiremen'.s  are  as 
follows: 

1  Processed  range  resolution,  Ar  (meters) 

2  Maximunr.  target  detection  range  (unambiguous),  (meters) 

3  Maximum  target  extent,  E  (meters) 

4  Target  radial  velocity,  v  (m/s) 

5  Allowable  range-doppler  coupling  e;  .jressed  in  terms  of  maximum  bin  sliifi  and 
dispersion 

1  he  adjustable  design  parameters  are  as  follows 

1  Frequency  step  size,  AF  (Hz) 

2  Pulse  width,  r  (sec) 

3  Number  of  pulses  writhin  coherent  processing  interv  al,  N 

4  Pulse  repehtion  frequency,  PRF  (Hz) 

5  Time  on  target,  Ti„.  (sec) 

The  design  process  consists  of  the  following  steps  in  the  order  listed  belovv 
1  From  the  required  processed  range  resolution,  calculate  the  waveform  bandwidth 

Kc 

=NAF  =  —  (3  24) 

2Ar 
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2  I'or  good  target  detectability,  determine  the  pulse  width  necessary  to  encompass  tl.e 
target 


c 

This  requiremen*  need  not  Oe  met  if  target  mapping  is  the  primary  concern 
J  Determine  the  fraction  of  the  unambiguous  range  window  which  is  occupied  by  the 
returns  from  an  onginal  range  bin  of  width  ct/2  This  fraction  is  giver  by  equation 
(3  15)  as  being  equal  to  the  product  tAF.  To  avoid  aliasing  of  the  profile,  tAF  cannot 
exceed  unity  In  the  absence  of  sfong  clutter,  this  product  should  be  set  at  o'"  near 
unity  to  provide  fine  processed  range  resolution 
4  The  frequency  step  size.  AF,  can  be  determined  by  simply  dividing  the  value  of  xAF 
by  the  pulse  width 


1 


(3.26) 


Provided  that  tAF  is  less  than  or  equal  to  unity,  the  frequency  step  size 
determined  by  equation  (3  26)  automatically  satisfies  the  criterion  that 


AF< 


c 

ZE 


(3  27) 


which  is  the  required  sampling  rate  in  the  frequency  domain 
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5  Determine  the  number  of  frequency  steps  N  based  on  the  required  bandwidth  and 


the  computed  frequency  step  size 


AF 


(3  28) 


6  Choose  the  PRF  as  follows 

A  For  stationary  targets,  select  the  PRF  based  on  the  following 
considerations 

i  Maximum  unambiguous  target  detection  range 


R 


uAvnb 


(3.29) 


ii  Transmitter  average  power 


~  ^(  fpRF  ) 


(3  30) 


B  For  moving  targets: 

i  Select  the  PRF  based  on  the  maximum  expected  target  radial  velocity 
and  the  maximum  tolerable  bin  shift, 


f  = 

*PRF 


f  vN 


B,„Ar(L_) 


(3  31) 
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2  Alternatively,  select  the  PRF  based  on  the  maximum  expected  target 
velocit)'  and  the  maximum  tolerable  dispersion  factor,  : 


f  = 
*PRF 


vN 


Ar(P„,„) 


(3.32) 


The  PRF  chosen  by  either  of  these  methods  will  be  consistent  since 
the  maximum  bin  shift  is  proportional  to  the  maximum  dispersion 
factor 


L 


max 


f 


(3  33) 


7  Determine  the  time  on  target  based  on  the  number  of  pulses  and  the  PRF 


(3  34) 


Upon  completion  of  these  steps,  one  must  check  to  ensure  that  the  chosen  PRF 
will  not  yield  an  unambiguous  range  in  conventional  radar  terms  (equation  3.29)  that  is 
les.s  than  the  maximum  detection  range  If  this  is  the  case,  then  the  PRF  must  be 
recomputed  in  step  five  with  a  different  value  of  maximum  target  radial  velocity  or  a  less 
stringent  maximum  bin  shift  factor  or  a  less  stringent  dispersion  factor 


3.  Design  Method  No.  3:  noise  environment,  fixed  time  oii  target 

This  design  method  is  based  on  the  assumptions  that  the  effects  of  clutt.  are 
negligible  and  the  time  on  target  is  fixed  by  other  design  considerations  such  a.  le 
antenna  beam  width  and  antenna  scan  rate  Moreover,  it  is  assumed  that  velocity 
compensation  has  not  been  applied  to  the  range  profile  The  principal  requirements  are  as 
follows: 

1  Processed  range  resolution,  Ar  (meters) 

2  Maximum  target  detection  range  (unambiguous),  (meters) 

3  Maximum  target  extent,  E  (meters) 

4  Time  on  target,  T,^  (sec) 

1  he  adjustable  design  parameters  are  as  follows 

1  Frequency  step  size,  AI'  (Hz) 

2  Pulse  width,  x  (sec) 

3  Number  of  pulses  within  coherent  processing  interval,  N 

4  Pulse  Repetition  Frequency,  PRF  (Hz) 

The  design  process  consists  of  the  following  steps  in  the  order  listed  below 

1  From  the  required  maximum  unambiguous  target,  calculate  the  PRF: 

(3  35) 

2  Calculate  the  numbei  of  pulses  N.  based  on  the  lime  on  target,  and  the  PRF 

N  =  T.,,(f,^-)  (3  36) 
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3  For  good  detectability,  detennine  the  pulse  width  necessary  to  encompass  the  target: 


(3.37) 


As  in  the  previous  design  method,  the  fraction  of  the  unambiguous  range  window 
occupied  by  the  original  range  cel!  should  be  at  or  near  unity,  in  the  absence  of 
strong  clutter.  However,  the  product  tAF  should  not  exceed  unity  to  avoid  aliasing 
4  Choose  the  frequency  step  size  AF  to  obtain  the  required  range  resolution  based  on 
the  number  of  pulses  N  determined  in  step  one 


AF  = 


Kc 

2NAr 


(3  38) 


The  major  difference  in  this  design  method  compared  to  the  previous  method  is 
that  the  PRF  is  determined  primarily  by  the  time  on  target  constraints  which  are 
determined  by  the  antenna  beam  width  and  scan  rate  As  a  result,  the  target  radial  velocity 
and  maximum  allowable  bin  shift  no  longer  play  a  role  in  determining  the  PRF  Moreover, 
the  number  of  pulses  within  the  burst  are  determined  by  unambiguous  range  and  time  on 
target  constraints  The  frequency  step  size  is  then  chosen  based  on  the  required  range 
resolution  and  the  precomputed  number  of  frequency  steps 

To  summarize,  the  primary'  design  parameters  of  a  step  frequency  radar  are  the 
number  of  pulses.  N,  the  frequency  step  size,  AF,  the  PRF,  and  the  pulse  width,  x 
Optima!  combinations  of  these  parameters  are  determined  to  achieve  certain  specified 
design  goals  Design  methods  are  developed  for  a  low  PRF  mode  of  operation  in  clutter 
limited  and  noise  limited  environments 
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IV.  SIMULATION  OF  HRR  PROFILES 


A  simulation  program  was  developed  using  MATLAB  to  graphically  show  the  effects 
of  various  step  frequency  waveform  parameters  on  a  synthetic  range  profile  The 
simulation  output  is  the  HRR  profile  corresponding  to  a  given  original  range  bin  of  width 
ct/2  There  are  three  versions  of  simulation;  one  for  noise  limited  situations,  one  for 
clutter  limited  situations,  and  one  for  clutter  limited  situations  which  includes  clutter 
cancellation  and  velocity  compensation  In  the  first  two  versions,  a  maximum  of  three 
moving  targets  within  the  original  range  cell  can  be  specified  The  third  version  simulates  a 
single  moving  target  Moreover,  all  major  waveform  parameters  such  as  the  number  of 
pulses  N,  the  frequency  increment  AF,  the  PRF,  and  the  carrier  frequency  can  be  specified 
In  addition,  the  signal  to  noise  ratio  (SNR)  and  the  signal  to  clutter  ratio  (SCR)  at  the 
output  of  the  IF  filter  can  be  specified  for  each  taiget  based  on  a  reference  value  for  a  1 
m’  target 

The  noise  source  is  modeled  as  white  gaussian  noise  with  a  variance  of  0  0002  The 
clutter  sources  are  modeled  as  a  series  of  ten  closely  spaced  sinusoidal  scatterers  with 
fixed  ranges  within  the  original  range  cell  The  magnitude  of  each  cluttei  source  is 
identical  and  is  based  on  the  signal  to  clutter  ratio  for  a  1  m*  target  A  complete 
description  of  all  three  versions  of  simulation  is  given  in  Appendix  A  The  program  codes 
are  listed  in  Appendix  B 
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A.  HRR  PROFILES  W  NOISE  LIMITED  SITUATIONS 


In  this  section,  the  HRR  profiles  will  be  embedded  in  white  gaussian  noise  with 
variai.ce  of  0.0002.  The  reference  waveform  parameters  are  listed  in  Table  4.1.  The 
reference  HRR  profile  is  shown  in  Figure  4  1 


TABLE  4  1  REFERENCE  PAR.AMETERS  FOR  HRR  SIMULATION  WITH  NOISE 


Figure  4  1  clearly  reveals  the  peaks  corresponding  to  the  three  targets  located  at  1204 
m,  1210  m,  and  1212  m  The  inverse  FFT  (IFFT)  bin  number  of  each  target  peak  is 
directly  related  to  the  target  range  To  determine  the  IFFT  bin  number  for  a  specific  target 
range,  the  foiiov.'ing  steps  are  earned  out  First,  the  modulo  range,  is  determined  by 
dividing  the  target  range,  R,,  by  the  unambiguous  range,  R^.,  as  follows 

Rmod  =  R.  "lod  R„  (4  1) 

The  EFFT  bin  number  corresponding  to  the  modulo  range  is  then  given  by 


IFFT,.,=  R„„, 


(4  2) 


In  this  case,  equals  15  meters  The  HRR  profile  in  Figure  4  1  represents  range  from 
1200  to  1215  meters 
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The  product  tAF  is  set  to  unity  in  the  reference  profile  which  results  in  the  original 
range  bin  of  width  ct/2  occupying  the  entire  unambiguous  range  window  .  This  gives  a 
processed  range  resolution  of  0.03  meters  unwindowed  However,  the  reference  pi  ofile 
has  been  multiplied  by  a  Hamming  window  to  reduce  the  sidelobes  of  each  target  peak 
which  results  in  a  slight  loss  of  range  resolution  from  0  03  m  to  0  04  m  In  the  following 
subsections,  the  pulse  vhdth,  t,  and  the  SNR  will  be  held  constant,  and  unless  specified 
other\\’ise,  all  HRR  profiles  are  multiplied  by  the  Hamming  window  The  remaining 
parameters  will  be  varied  one  at  a  time  from  their  reference  value,  and  the  resulting  effects 
on  the  range  profile  will  be  displayed 

i.  Effect  of  Varying  the  Number  of  Pulses 

Figure  4  2a  shows  the  effect  of  increasing  the  number  of  pulses  to  1000 


The  major  change  in  increasing  the  number  of  pulses  is  the  increased  range  resolution 
which  is  now  0  015  m  (unweighted  )  Moreover,  there  is  an  increase  in  the  SNR  due  to  the 
increase  in  processing  gain  due  to  coherently  integrating  1000  instead  of  500  pulses  It  the 
number  of  pulses  is  reduced  from  the  original  number,  the  range  resolution  becomes 
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coarser,  and  the  SNR  decreases  due  to  a  reduction  in  processing  gain  These  results  are 
shown  in  Figure  4.2b.  In  this  figure,  the  number  of  pulses  has  been  reduced  to  250  wiiich 
results  in  a  range  resolution  of  0.06  m  The  reduction  in  the  number  of  pulses  reduces  the 
SNR  gain  to  such  an  extent  that  the  third  scatterer  is  not  detectable  from  noise 


Figure  4  2b  HRR  Profile  of  three  stationary  targets  in  white  noise  N  =  250 


Note  the  third  peak  (bin  200)  is  not  distinguishable  from  the  noise 

2.  Efl'ect  cf  a  Shift  in  Target  Range 

One  of  the  key  features  apparent  from  the  HRR  profile  is  that  it  is  highly 
ambiguous  in  absolute  range  Only  relative  scatterer  position  within  the  original  pulse 
width  is  indicated  in  the  profile  Targets  located  at  multiples  of  the  unambiguous  range 
would  all  lie  in  the  same  relative  position  in  the  HRR  profile  In  addition,  a  slight  shift  in 
range  position  of  the  targets  relative  to  the  radar  will  cause  the  peaks  to  shift,  and  possibly 
wraparound,  resulting  in  range  ambiguities  For  example.  Figure  4  3  shows  the  range 
profile  when  the  three  targets  each  have  shifted  outward  by  six  meters  in  range  As  a 
result,  the  second  and  third  peaks  in  the  original  profile  have  wrapped  around  so  that  they 
appear  first  in  the  profile,  whereas  the  first  peak  in  the  original  profile  appears  last 
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3.  Effect  of  Varying  the  Frequency  Step  Siit 

The  frequency  step  size  determines  not  only  the  range  resolution,  but  also  the 
size  of  the  unambiguous  range  v.-indov.’,  A  larger  value  of  AF  results  in  increased  range 


resolution  and  a  smaller  unambiguous  range  window.  Figure  4  4a  shows  the  effect  of 
increasing  the  step  size  to  20  MHz  The  unambiguous  range  has  been  reduced  to  7  5 
meters,  and  the  range  resolution  was  increased  to  0.015  m  Therefore,  the  peaks  at  initial 
ranges  ot  1210  m,  and  1212  m  have  been  folded  over  The  third  peak  is  not  resolvable 
because  its  relative  location  is  now  very  close  to  the  sidelobes  of  the  largest  peak  This 
figure  indicates  a  situation  in  which  the  unambiguous  range  window  is  smaller  than  the 
original  range  bin  width  cf  ct/2,  due  to  undersampling  in  the  frequency  domain 

Figure  4  4b  shows  the  opposite  situation,  in  which  the  frequency  step  size  has 
been  reduced  to  5  MHz  In  this  case,  the  range  resolution  decreases  to  0  06  m,  but  the 
width  of  the  unambiguous  range  window  increases  to  30  m  This  figure  indicates 
oversampling  in  the  frequency  domain,  because  the  entire  contents  of  one  original  range 
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cell  of  width  ct/2  are  contained  in  the  first  half  of  the  window  The  second  half  represents 
a  free  region  clear  of  any  targets  and  possible  clutter,  and  contains  only  background  noise 


I 


Figure  4  4b  HRR  Profile  of  thiee  stationary  targets  in  white  noise  with  AF  =  5  MHz 


=  30  m  The  entire  contents  of  one  original  range  bin  now  comprise  one 
half  of  the  unambiguous  range  window 
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4.  Effect  of  Varying  the  PRF  and  the  Carrier  Frequency 

For  stationary  targets,  the  pulse  repetition  frequency  (PRF)  has  no  effect  on  the 
final  FIRR  profile  However,  dispersion  and  bin  shifting  of  profiles  are  sensitive  to  PRF, 
and  these  only  occur  if  the  target  has  non-zero  radial  velocity.  In  addition,  changing  the 
nominal  carrier  frequency,  f^,  has  no  effect  on  the  range  profile  of  stationary  scatterers 
Figure  4.5  shows  the  HRR  profile  for  three  stationary  targets  when  the  PRF  is  increased 
to  200  kHz,  and  Figure  4.6  shows  the  range  profile  for  a  carrier  frequency  of  12  GHz 


Figure  4  5  HRR  Profile  of  three  stationary'  targets  in  white  noise  PRF  =  200  kHz 


This  figure  is  identical  with  Figure  4  1  which  displays  the  results  with  a 
PRF  of  50  kHz 


frequency  set  to  12  GHz  The  results  are  nearly  identical  with  Figure  (4  1) 
which  showed  the  results  for  a  carrier  frequency  of  1  GHz 

5.  Effect  of  Range  Weighting  Function 

Generation  of  HRR  profiles  for  a  given  scatterer  via  the  step  frequency 
w'aveform  will  result  in  range  sidelobes  which  can  degrade  the  ability  to  detect  a  target  in 
close  proximity  to  another  target  The  inverse  FFT  results  in  a  5in(x)/x  shape  for  each 
scatterer  peak  which  results  in  a  peak  sidelobe  level  of -13  dB  (-43  dB  for  Hamming  [Ref 
3:p  426]  )  below  the  main  lobe  The  application  of  a  range  weighting  function  prior  to 
taking  the  inverse  FFT  will  result  ir;  reduced  peak  sidelobe  level  The  reduction  of 
sidelobe  levels  comes  at  the  expense  of  a  slight  reduction  in  the  peak  signal  level  and  a 
reduction  in  range  resolution  due  to  a  broadening  of  the  main  peak  [Ref  6:pp.  246-249] 
The  reduction  in  resolution  is  characterized  by  a  dimensionless  parameter,  K,  v/hich 
appears  in  the  expression  for  the  processed  range  resolution,  equation  (2  3)  Figure  4  7a 
shows  an  HRR  range  profile  without  a  weighting  function  whereas  Figure  4  7b  shows  the 
same  profile  multiplied  by  a  Hamming  window  The  range  resolution  for  the  unweighted 
profile  is  0  03  m  as  compared  to  0  039  m  for  the  windowed  profile 
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Magnitu<l«  response  in  dB 


6.  EfTect  of  Target  Vdocity 


Target  motion  degrades  the  resolution  of  the  HRR  profile  in  three  ways  First, 
the  target  peak  is  broadened  into  several  range  bins.  The  number  of  processed  range  bins 
that  the  dispersed  peak  occupies  is  equal  to  the  dispersion  factor  P  given  by 


P  = 


Nv 

fpRf^ 


(43) 


Moreover,  the  dispersion  of  a  target  peak  will  also  result  in  a  reduction  in  SNR  An 
extremely  broad  target  peak  will  also  reduce  the  detectability  of  other  close-in  targets 
Finally,  target  motion  will  lead  to  a  shift  in  the  target  peak  location  The  number  of  bins 
sl'iifted  is  equal  to  the  bin  shift  factor  L  given  by 


L  = 


^  f 

*c 


J 


Nv 


(44) 


If  the  bin  shift  factor  L  exceeds  the  diffen^nce  between  the  number  of  pulses  N  and  the 
IFFT  bin  location  for  a  corresponding  stationary  target,  the  target  peak  will  wrap  around 
resulting  in  an  aliased  profile 

The  effects  of  target  i..v  iicn  can  be  mitigated  by  either  increasing  the  PRF  or  by- 
velocity  compensation  An  increase  in  the  PRF  will  decrease  the  dispersion  factor  P  and 
the  bill  shift  factor  L  as  in  equations  (4  3)  and  (4  4)  respectively  Switching  to  a  higher 
PRF  will  improve  the  range  resolution  and  detectability  of  moving  targets  without 
requiring  prior  knowledge  of  the  target  velocity  However,  a  higher  PRF  results  in  a 
reducer  ambiguous  range  (in  conventional  radar  terms) 


2f 


PRF 


(4.5) 


Velocity  compensation,  on  the  other  hand,  converts  the  rarge  profile  of  a  moving  target 
into  the  corresponding  profile  of  a  stationary  target  at  the  initial  range  location  As  a 
result,  the  dispersion  and  bin  shifting  are  completely  canceled  However,  velocity 
compensation  requires  an  accurate  estimate  of  the  target  velocity  If  different  scatterers 
within  one  pulse  width  have  different  velocities,  an  average  of  the  velocities  must  be  used, 
but  it  will  not  result  in  complete  compensation  fcr  any  scatterer. 
a.  Effect  of  Uncompensated  Target  Velocity 

Figures  4.8a-d  show  the  uncompensated  HRR  profile  for  two  stationary 
targets  and  one  moving  target  with  speeds  of  100  m/'s,  300  m/s,  500  m/s,  and  -200  m/s 
respectively,  at  a  PRF  of  50  kHz 
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Figure  4.8b  HRR  Profile  of  one  moving  target,  (T1 )  wuh  radial  velocity  of  300  m/sec 

PRF  =  50  kHz  T2  and  T3  are  Stationary- 


Figure  4.8c  HRR  Profile  of  one  moving  target.  (T1 )  with  radial  velocity  of  500  'n/sec 


PRF  =  50  IcFIz  T2  and  T3  are  stationary 


50 


T2  and  T3  are  stationai^ 

Figure  4  9  shows  the  effect  of  all  three  scatterers  with  different  velocities 
with  respect  to  the  radar.  This  situation  can  exist  if  different  regions  of  a  target  have 
differentia!  velocities  with  respect  to  each  other.  Figure  4. 10  shows  the  range  profile  it  all 
three  scatterers  have  the  same  radial  velocity.  The  radar  cross  section  of  each  moving 
target  has  been  set  to  10  m’  in  Figures  4.1 1  to  4  13  for  better  visual  clanty  in  identifying 
the  effects  of  multiple  moving  targets  on  the  range  profile 


5! 


Figure  4  9  HRR  Profile  for  three  moving  targets  each  with  a  different  radial  velocity 


The  radial  velocities  are  1 50  m/s  (Tl),  100  m/s  (T2),  and  50  m/s  (T3) 
PRF  =  50  kHz 


Figure  4  10  FIRR  Profile  of  three  targets  each  with  a  radial  velocity  of  100  m./sec 


PPJ  =  50  kHz 
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h.  Effect  of  Velocity  Compensation 


Velocity  compensation  at  a  speed  of  100  m/sec  can  be  applied  to  the  profile 
in  Figure  4. 1C  to  enhance  its  range  resolution  whereas  the  same  compensation  factor  will 
not  achieve  optimum  results  for  Figure  4  9  which  has  different  radial  velocities.  Figures 
4  1 1  and  4.12  display  the  results  of  velocity  compensation  applied  to  Figures  4.9  and  4.10 
respectively 


0  WS  255  iCb  'nsJo 

IFrr  b»n  numbei 


Figure  4  1 1  HRR  Profile  in  Figure  (4.9)  with  velocity  compensation  of  1 00  nVs 

Note  that  T2  (at  100  m/s)  is  compensated  correctly,  but  that  Tl,  and  T3 
are  still  dispersed  due  to  a  mismatch  between  the  target  velocity  and  the 
compensated  velocity 
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c.  Effect  of  Increasing  the  PRF 

The  major  limitation  of  applying  velocity  compensation  factors  is  that  an 
accuiate  estimate  of  the  velocities  of  individual  scatterers  must  be  obtained  a  priori 
Moreover,  differential  velocities  between  scatterers  within  a  single  target  wi'!  result  in 
imperfect  compensation  if  a  single  velocity  correction  is  applied,  as  shown  in  Figure  4  11 
As  a  result,  the  PRF  luould  be  increased  as  much  as  possible  within  other  design 
constraints  to  reduce  the  negative  impact  of  radial  velocity  on  the  range  profile  As  the 
PRF  is  increased,  the  velocity  correction  becomes  less  stringent  Figures  4  13a-c  show  the 
effect  of  varying  the  PRF  from  20  kHz  to  100  kHz,  and  250  kHz  respectively  for  three 
scatterers  each  having  a  radial  velocity  of  100  m/'s 
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Figure  4  13a  HRR  Profile  of  three  moving  targets  with  a  radial  velocity  of  100  m/sec 


PRF  =  20  kHz.  Note  the  increased  dispersion  and  bin  shift  of  all  targets, 
and  the  aliasing  caused  by  the  wraparound  of  the  target  peak  (T3) 


•100 _ _ . _ _ _ 

0  100”  ~^00  500  400  500 
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Figure  4  13b  HRR  Profile  of  three  targets  with  a  radial  velocity  of  100  nv'sec 
PRF  =  100  kHz  Note  the  reduction  in  dispersion  and  bin  shifting 
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Figure  4  13c  HRR  Profile  of  three  scatterers  with  radial  velocity  of  1 00  m/s 


PRF  =  250  kHz.  Note  that  this  profile  looks  similar  to  the  velocity 
compensated  profile 

B.  HRR  PROFILES  IN  CLUTTER  LIMITED  SITUATIONS 

!f  echo  returns  from  clutter  are  present  in  the  original  range  cell,  then  the  appeal auce 
of  the  HRR  profile  changes  dramatically.  The  processed  range  p.’-ofilc  contains  numerous 
peaks  due  tc  clutter  as  well  as  a  few  peaks  due  to  the  targets  of  interest 

By  keeping  the  product  of  xAI-  <  > ,  the  clutter  filled  range  cell  occupies  only  a  small 
percentage  of  the  unambiguous  range  window  depending  upon  the  value  of  xZiF  as  given 
by: 

=  tAF  (4  5) 

The  remainder  of  the  range  window  consists  of  a  cluiter-free  legion  containing  only 
thermal  noise  Targets  having  sufficient  radial  velocity  in  the  line  of  sight  of  the  radar  will 
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have  their  peaks  shifted  from  the  clutter  filled  region  to  the  clutter-free  region  of  the  range 
profile  where  it  wili  be  easy  to  detect  them  The  extent  to  which  the  moving  target 
migrates  into  the  clutter-free  region  depends  on  the  parameters  which  determine  the 
dispersion  factor  P  and  the  bin  shift  factor  L  as  given  by  equations  (4.3)  and  (4  4) 

Application  of  velocity  compensation  in  a  clutter  limited  situation  without  clutter 
cancellation  is  counterproductive  since  the  compensation  factor  causes  the  target  peak  to 
move  back  into  the  clutter  region.  Although  the  dispersion  of  the  target  peak  results  in  a 
loss  of  resolution,  the  moving  target  will  be  detected  whereas  the  corresponding  stationary' 
target  must  compete  with  the  clutter  return  in  order  to  be  detected 

The  HRR  profile  of  three  scatterers  within  a  single  original  range  cell  in  the  presence 
of  clutter  is  simulated  by  adding  several  stationary'  peaks  to  represent  the  clutter.  The 
nominal  value  of  tAF  is  chosen  to  be  0. 1  Therefore,  the  clutter  region  of  the  range  profile 
comprises  ten  percent  of  the  unambiguous  range  window.  The  signal  to  clutter  ratio 
(SCR)  is  the  ratio  of  the  signal  power  to  the  clutter  power  at  the  input  of  the  signal 
processing  unit  for  a  1  m*  target  The  clutter  magnitude  is  also  based  on  the  amplitude  of 
a  1  m'  target  The  SCR  can  be  specified  independent  of  the  SNR  by  the  user  The  full 
description  of  the  simulation  program  is  listed  in  Appendix  A.  The  parameters  for  the 
reference  HRR  profile  with  clutter  are  shown  in  Table  4  2  The  corresponding  profile  is 
displayed  in  Figure  4  14  Unless  mentioned  otherwise,  all  other  parameters  except  the 
parametei  of  interest  will  have  the  same  value  as  shown  in  the  reference  profile 
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TABLE  4  2  REFERENCE  PARAMETERS  FOR  HRR  PROFILE  WITH  CLUTTER 


Parameter 


Number  ol  Pulses  (N) _ 

Frequency  Step  Size  (AT) 


Pulse  Repetition  Frequency  (PRF) 
Pulse  Width  (t 


Nominal  Carrier  Frequency  (f) 


Unambiguous  Range  Window  (P  J 


Range  Weighting  Function _ _ 

Initial  Target  Ranees 


Reference  Value 


500 


I  MHz 


20  kHz _ 

0  1  usee 


1  GHz 


150  in 


Hamniin 


reel 

SC 

R  ratio 

.  -20dB(T3 

,  lU  Ott 


Target  Velociti-'s 


'■1 

•20  T3 

I  I  1 


-^0 1 1  L  I. 


S.  -60 
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Figure  4  !4  Re!erence  HRR  Profile  for  three  s'.at'onary  targets  in  clutte. 


Due  to  clutter  in  tl,^  original  range  bin,  only  the  two  large,  stat’onary  targets  are 
above  the  clutter.  The  clutter  region  comprises  a  total  of  (tAF)N  bins,  which  in  this  case 


equals  50. 

1 .  Effect  of  Varying  the  Number  of  Pulses 

An  increase  in  the  number  of  pulses,  N,  results  in  increased  range  resolution 
compared  to  the  resolution  of  0  .3  m  (unwindowed)  for  Figure  4,14  The  effects  of 
increasing  the  number  of  pulses  to  1000  and  reducing  the  number  of  pulses  to  250  are 
shown  in  Figure  4. 1  .Sa,  and  Figure  4  1 5b  respectively 
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-100 


0  ^05  ^  650^^ - 

!FFT  bin  number 

■igure  4  1 5a  HRR  Profile  for  three  stationary’  targets  in  clutter  N  =  1000 
Note  the  increased  range  resolution  of  0  1 5  m  (unwindowed) 


Figure  4. 15b  HRR  Profile  for  three  stationary  targets  in  clutter.  N  =  250 
Note  the  reduction  in  range  lesolutioii  (0  6  m  unwindowed) 

2.  Effect  of  Varying  the  Frequency  Step  Size 

An  increase  in  the  frequency  step  size  also  results  in  finer  range  resolution 
Howevei  ,  the  clutter  occupies  a  larger  fraction  of  the  unambiguous  range  window  as 
shown  in  Figure  4  16.  As  a  result,  there  is  less  clutter  free  space  with  which  to  distinguish 
between  moving  targets 
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resolution  =  0  06  m  (unwindowed).  Note  that  the  clultei  occupies  half  of 
the  unambiguous  '»r.ge  window 


3.  Ef.ect  of  Target  Velocity 

In  a  clutter  limited  situation,  targets  with  sufficient  velocity  will  rrJgrate  into  the 
clutter  free  legion  Assuming  that  the  number  of  pulses  N  is  held  fixed  to  maintain  the 
required  range  resolution,  the  PRF  primarily  determines  the  extent  to  which  the  target 
moves  into  the  clutter  free  region  As  given  by  equations  (4.3)  and  (4.4),  the  dispersion 
factor  and  bin  shift  factor  are  inversely  proportional  to  the  PRP  Thu.s  for  a  given  target 
velocity,  a  lower  PRF  will  shift  the  target  peak  further  into  the  clutter  free  region 
a.  Effect  of  Single  Moving  Target 

The  eftects  of  velocity  changes  for  the  100  m'  target  (Tl)  are  displayed  in 
Figures  4  l/a-d  The  increase  in  radial  velocity  from  300  m/s  to  1000  m/s  results  in  both 
increased  dispersion  of  the  scatierer  peak  as  well  as  an  increased  amount  of  bin  shifting  of 
the  profile  An  incoming  target,  as  shown  in  Figure  4  17d  shows  the  same  effects  except 
that  the  bin  shift  will  be  in  the  opposite  direction 
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b.  Effect  of  Multiple  Moving  Targets 

The  effects  of  having  multiple  moving  targets  with  the  same  velocity  are 
shown  in  Figures  4  18a  and  (4.18b)  Tliis  situation  could  correspond  to  multiple  scatterers 
within  one  distributed  target  traveling  at  the  same  velocity.  As  the  velocity  is  increased 
from  100  m/s  to  500  m/s,  the  dispersed  target  peaks  merge  into  each  other  making 
individual  target  detection  nearly  impossible 
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'igure  4  1 8a  HRR  Profile  of  three  moving  targets  in  clutter  Radial  velocity  =  !  00  nx^s 
Note  that  each  dispersed,  shifted  target  peak  is  distinguishable 
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Figure  4  18b  HRR  Profile  of  three  moving  targets  in  clutter  Radial  velocity  =  500  m/s. 
Note  the  merging  of  the  target  peaks  into  one  large  peak 


The  effects  of  multiple  moving  targets  with  different  velocities  on  the  range 
profile  are  presented  in  Figures  4  19a,  and  (4  19b)  This  situation  can  indicate  relative 
motion  between  diffeieru  scattering  centers  within  a  larger  distributed  target  or  closely 
spaced  distinct  targets  As  shown  in  Figure  4  19b,  it  is  possible  to  mistake  a  ",si  opening 
target  with  a  rapidly  closing  target  Therefore,  moving  targets  in  the  left  half  of  the  profile 
should  be  interpreted  as  opening  targets  (negative  Doppler),  and  targets  in  the  right  haif 
should  be  viewed  as  closing  targets  (positive  Doppler )  In  both  cases,  the  magnitude  of 
the  Doppler  shift  increases  as  the  targets  approach  the  center  The  sign  convention  is  the 
opposite  of  that  used  in  conventional  pulse-doppler  radars 
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Figure  4  1 9a  HRR  Profile  of  three  moving  targets  in  clutter  Radial  velocities  for  T 1 , 


T2,  and  T3  are  150  m/s,  300  m/s  and  450  m/s  respectively 


Figure  4  19b  HRR  Profile  of  three  moving  targets  in  clutter  T1  is  an  incoming  target 


with  a  radial  velocity  of  -500  nx's.  T2  has  a  radial  velocity  of  450  nvs,  and 
T3  has  a  radial  velocity  of  100  m's  Since  T1  and  T2  have  nearly  opposite 
velocities,  it  is  possible  to  mistake  a  fast  opening  target  (T2),  for  a  fast 
closing  target  (Tl) 


c.  Effect  of  Varying  the  P.\F 


The  effects  of  var>ing  the  radar's  pulse  repetition  frequency  are  shown  in 
Figures  4  20a-c  for  PRF's  of  5  kHz,  20  kHc.  and  200  kHz  respectively  for  three  targets 
with  velocities  of  150  m/s,  300  m/s,  and  450  m/s  For  these  plots,  the  cross  sections  of  all 
targets  have  been  set  at  10  m^  to  focus  solely  on  the  effects  of  varying  the  PRF  For  the 
low  PRF  case,  the  bin  shift  and  dispersion  of  the  third  target  result  in  fold  over  of  the 
peak,  resulting  in  an  aliased  profile  As  the  PRF  is  increased  to  200  kHz,  the  peaks 
broaden  slightly  but  remain  in  the  clutter  region 


Figure  4,20a  HRR  Profile  of  three  moving  targets  in  clutxer  PRF  =  5  kHz  Note  the 


increased  bin  shift  and  dispersion  of  the  target  peaks  and  the  foldover  of 
the  third  target 


Magnitude  response  in  d8 


Figure  4  20b  HRR  Profile  of  three  moving  targets  in  clutter.  PRF  =  20  kHz 


b  TB6  job  300  405  '500 
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Figure  4  20c  HRR  Profile  of  three  moving  targets  in  clutter  PRF  =  200  kHz 

Note  that  all  three  targets  are  still  in  the  clutter  due  to  the  high  PRF 
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Finally,  velocity  compensation  cannot  be  applied  to  moving  targets  if  clutter  is 
present  in  the  original  range  cell.  Applying  a  compensation  factof  will  eliminate  the 
dispersion  of  the  peak  resulting  in  increased  range  resolution.  However,  unless  the  clutter 
is  canceled,  the  compensation  will  shift  the  target  peak  back  into  the  clutter  Figure  4.21 
shows  the  effect  of  velocity  compensation  for  tliree  moving  targets  each  with  a  radial 
velocity  of  100  m/s  at  a  PRF  of  20  kHz  As  a  result  of  the  compensation,  each  of  the 
targets  have  shifted  back  into  the  clutter  region  and  are  less  detectable 
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figurt  4  21  HRR  Profile  of  three  moving  targets  in  clutter  Radial  velocit>  =  100  m/s 


for  each  target  Velocity  compensation  of  100  ni/s  has  been  applied  to  the 
profile  Note  that  the  targets  have  sh-  fed  back  i;Ko  the  clutter 
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C.  IMPROVING  RESOLUTION  OF  MOVING  TARGETS  BY  CLUTTER 
CANCELLATION  AND  VELOCITY  COMPENSATION 

The  range  profiles  of  moving  targets  in  clutter  are  clearly  visible  only  i^"  the  target  has 

sufficient  velocity  to  migrate  to  the  clu‘ter  free  region  -Although  the  target  peak  will  be 

easily  detected  once  in  the  clutter  fVce  rone,  the  resulting  resolution  v»il!  be  poorer  due  to 

dispersion  The  dispersion  factci  P  given  in  equation  (4  3)  determines  the  number  of 

processed  range  bins  occup'ed  by  the  moving  target  The  etfective  processed  resolution  of 

the  target  equals  PAr  instead  of  the  optimum  value  of  Ai  In  addition  to  the  loss  of 

resolution,  an  excessive  shirt  in  the  peak  location  due  to  range-doppler  coupling  may 

result  in  target  wrap  around  The  number  of  bins  shifted  is  equal  to  the  bin  shift  factor  L, 

given  in  equation  (4  4) 

To  overcome  these  drawbacks  one  could  apply  velocity  compensation  directly 
However,  this  would  only  shift  the  target  peak  back  into  the  clutter  Therefore,  the  clutter 
must  be  canceled  prior  to  applying  the  velocity  compensation  One  method  to  impro''e  the 
range  resolution  of  moving  targets  in  clutiei*  is  proposed  as  follows  First,  the  frequency 
domain  target  reflectivity  data  is  converted  to  the  time  domain  HRR  profile  by  an  inverse 
FFT  The  frequency  domain  data  is  not  windowed  prior  to  the  IFFT  The  clutter  filled 
region  of  the  range  profile  is  canceled  by  multiplication  with  a  rectangular  gating  function 
The  modified  range  profile  is  then  windowed  prior  to  being  transformed  into  frequency 
domain  data  via  the  FFT  The  velocity  compensation  factor  is  applied  to  the  frequency 
dom..in  data,  and  the  resulting  data  is  transformed  to  the  corrected  range  profile  by  an 
inverse  FFT  The  steps  are  summarized  below 

1  Convert  frequency  domain  reflectivity  data  to  the  initial  HRR  profile  via  the  ITFT 

2  Apply  clutter  cancellation  to  the  initial  HRR  profile 

3  Estimate  the  compensation  velocity 

4  Apply  a  windowing  function  to  the  initial  HRR  profile 
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5  Convert  the  modified  HRR  profile  to  frequency  domain  data  via  the  FFT 

6  Apply  velocity  compensation  to  the  resulting  frequency  domain  data 

7.  Convert  the  compensated  frequency  domain  data  to  the  final  HRR  profile  by  EFFT. 

As  a  result  of  applying  these  steps,  the  moving  target  peak  moves  to  the  correct  range  in 
the  HRR  profile  More  imponantly,  the  resulting  processed  range  resolution  equals  Ar 
•nstead  of  PAr 

One  key  assumption  in  this  method  is  that  all  targets  within  the  original  range  bin 
have  the  same  velocity  This  assumption  is  valid  in  the  case  of  several  distinct  scatterers 
origin '^ting  from  a  large  extended  target  As  a  result,  a  single  velocity  compensation  factor 
can  be  applied  to  the  HRR  profile  corresponding  to  an  original  range  bin  Since  the  target 
\'clocity  is  unknown,  a  number  of  different  velocity  compensations  can  be  applied  one  at  a 
time,  and  the  compensation  wliich  yields  the  best  resolution  can  be  taken  as  the  target 
veloc’ty 

A  simulation  program  was  developed  in  MATLAB  to  demonstrate  clutter 
cancellation  with  velocity  compensation  for  a  single  moving  target  This  program  requires 
the  user  to  supply  the  IFFT  bin  number  corresponding  to  the  maximum  extent  of  the 
flutter  region,  and  the  compensation  velocity  The  reference  parameters  for  the  simulation 
are  listed  in  Table  4  3  A  complete  description  of  the  program  (HRRCAD  M)  is  given  in 
Appendix  A  As  opposed  to  earlier  profiles,  the  range  profiles  are  given  in  magnitude 
instead  of  decibel  units  for  improved  visual  clarity  All  simulation  parameters  other  than 
that  wnich  is  being  varied  are  fixed  at  the  reference  values  indicated  in  Table  4  3 


TABLE  4  3  REFERENCE  PARAMETERS  FOR  HRR  SIMULATION  WITH 

CLUTTER  CANCELLATION  AND  V'ELOCITY  CON4PENSATION 


- - — ^ — 

Faraincler 

Reference  Value  ! 

Number  of  Pulses  (N) 

500 

Frequency  Step  Size  (AF) 

1  MHz 

Pulse  Repetition  Frequency  (^RF) 

20  kHz 

Pulse  Width  (t) 

0  1  psec 

1  Pr'^cessed  Rang"  Resoiutio"  (Ar) 

0.3  m 

Nominal  Carrier  Frequency  (f„) 

1  GHz 

Unambiguous  Range  Window  (RJ 

150  m 

Fraction  of  R„  occupied  by  original 
ranee  bin  (tAF) 

0  1 

' 

Range  Weighting  Function 

Hamming  (K=  1  33) 

1  _ 

111 

500  m/s 

1  Target  Cross  Sec  lon 

1  m' 

1 

1  SCR  for  1  m'  target 

.20®  _  -  i 

I.  Effert  of  Compensation  Velocit)' 

The  eflfect  of  applying  velocity  compensation  factors  of  0  ni/s  (uiico'^.^iensaied), 
100  n'Js,  300  rrv's,  300  ni's,  700  ni^s  and  1000  m/'s  for  the  reference  target  after  applying 
clutter  cancellation  is  shown  in  Figures  4  22a-f  respectively  In  each  of  these  profiics,  the 
effects  of  noise  are  assumed  to  be  negligible  compared  to  cluttei  All  other  reference 
parameters  are  held  constant 
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Magnitude  response 


Figure  4  22a  Uncompensated  HRR  Profile  Target  velocity  =  500  ni/s 
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0  100  200  300  ■  400  —500 

IFFT  810  No 

'igure  4  22b  HRR  Profile  with  clutter  cancellation  Target  velocity  =  500  m/s 
Compensation  velocity  =  1 00  m/s 
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Figure  4  22c  HRR  Profile  with  after  cancellation.  Target  velocity  =  500  ni/'s 
Compensation  velocity  =  300  ni/s. 


IFF1  Bin  No 


Figure  4  22d  lIRR  F/of.le  v>ith  cF'tter  cancellation  Target  velocity  =  500  ra/s 

Compensation  velocity  -•  500  nv's  (i  e  correct  velocity  compensation) 
Note  the  improvement  in  range  resolution  for  the  target  peak 
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CompanMitM  HRR  pnini* 

01* _ _ _ _ _ _ _ 


IFHT  Bin  No 

Figure  4  22e  FCRJR.  Profile  with  clutter  compensation  Target  velocity  =  500  m/s 
Compensated  velocity  =  700  m/s 
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Figure  4  22f  HRR  Profile  with  clutter  cancellation  Target  velocity  =  500  m  s 
Compensated  velocity  =  1000  m/s 


As  the  velocity  compensation  factor  approaches  the  tnie  target  velocity,  the 
target  peak  sharpens,  resulting  in  increased  signal  strength  and  finer  range  resolution.  This 
technique  will  increase  the  probability  of  detectioi,  as  the  output  is  higher  for  conect 
compensation  It  is  also  pr'ssible  to  determine  the  target  velocity  as  the  conect  velocity 
compensation  gives  a  larger  output  as  compared  with  incorrect  compensations 
2.  Effect  of  Compeosation  Velocity  at  a  Reduced  Signal  to  Clutter  Ratio 

The  effects  of  varys«ig  the  compensation  velocity  from  0  m/s  to  100  m/s,  300 
m/s,  500  m/s,  700  m/s,  and  1000  m/s  :t  an  SCR  of  -20  dB  are  displayed  in  Figures  4  23a-f 
respectively  From  these  figutes,  it  is  evident  that  the  signal  strength  is  maximized  when 
the  compensation  velocity  equal?  the  actual  taiget  vebcitv  These  HRR  profiles  fo'low  the 
same  pattern  as  those  vdth  a  high  SCR  of  20  dB  which  were  displayed  in  Figures  4  22a-f 
The  major  difference  is  in  the  relative  amplitudes  of  the  target  peaks 


Figure  4  23b  HRR  Profile  with  clutter  cancellation  SCR  =  -20  dB 


Target  velocity  =  500  m/s  Compensated  velocity  =  100  m/s 


Figure  4.23c  HRR  Profile  with  clutter  cancellation  SCR  =  -20dB 


Target  velocity  =  500  m/s  Compensated  velocity  =  300  nVs 
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Figure  4.23d  HRR  Profile  wth  cluner  cariceliation.  SCR=^-20dB 

Target  velocity  =  500  m/s  Compensated  velocity  =  500  m/s 


I 


Figure  4  23e  HRR  Profile  with  clutter  cancellation  SCR='-20dB 

Target  velocity  =  500  m/s  Compensation  velocity  700  ni/s 
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Figure  4  23f  HRR  Profile  with  clutter  cancellation.  SCR  =  -20  dB. 

Target  velocity  =  500  m/s  Compensation  velocity  =  1000  m/s 

In  summary,  a  simulation  program  was  developed  to  display  the  HRR  profile 
corresponding  to  a  single  original  range  bin  of  width  ct/2  There  are  three  versions  of  this 
program,  one  for  noise  limited  environments,  one  for  clutter  limited  environments,  and  one 
for  clatter  limited  environments  wliich  includes  clutter  cancellation  followed  by  repeated 
velocity  compensations  The  effects  of  individually  varying  waveform  parameters  such  as 
the  number  of  pulses,  N,  the  frequency  step  size,  AF,  and  the  PRF  on  the  HRR  profile  are 
graphically  illustrated,  as  well  as  the  effects  of  target  motion  and  velocity  compensation 


V.  CONCLUSIONS 


This  thesis  investigates  the  realization  of  high  radar  range  resolution  by  means  of  the 
step  frequency  waveform  This  waveform  achieves  high  range  resolution  by  coherently 
processing  the  returns  from  N  pulses  each  having  a  different  carrier  frequency  that 
changes  by  a  fixed  amount,  from  pulse  to  pulse  The  step  frequency  waveform  has  two 
principal  advantages  compared  to  conventional  high  resolution  waveforms.  First,  since  the 
effective  waveform  bandwidth  is  realized  over  several  pulses,  each  individual  pulse  is 
sufficiently  wide  which  eases  the  A/D  sampling  requirements  Second,  the  step  frequency 
waveform  can  be  implemented  by  conventional  radar  equipment,  with  the  addition  of  a 
stepped  frequency  synthesizer  Thus  high  range  resolution  can  be  obtained  using 
narrowband  components.  The  primary  limitation  of  the  step  frequency  waveform,  as 
Compared  to  intrapulse  pulse  compression,  is  the  time  required  to  transmit  and  receive  N 
pulses  needed  to  obtain  the  equivalent  v^de  bandwidth  o  a  single  narrow  pulse 

The  effects  of  various  waveform  parameters  on  the  high  resolution  range  profile 
(HRR  profile)  are  determined  for  both  stationary  and  moving  targets  in  noise  and  clutter 
environments  The  primary  parameters  include  the  number  of  pulses  N,  the  frequency  step 
size  AF,  the  PRF,  and  the  pulse  width  t  Design  methods  are  developed  based  on  specified 
requirements  to  optimize  system  performance  for  a  low  PRF  mode  in  clutter  limited  and 
noise  limited  situations  The  design  process  specifies  the  range  of  values  of  each  of  the 
primary  waveform  parameters  which  results  in  the  required  system  performance 

Target  motion  in  the  presence  of  noise  or  clutter  results  in  dispersion  of  the  target 
peak  in  the  HRR  profile  resulting  in  lower  SNR  as  well  as  a  shift  in  the  peak  location  due 
to  range-doppler  coupling,  which  can  also  result  in  target  wraparound  To  overcome  these 
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effects,  velocity  compensation  can  be  applied  to  the  range  profile  provided  that  an 
accurate  target  velocity  estimate  is  available  If  velocity  compensation  is  not  possible, 
increasing  the  PRF  can  mitigate  the  effects  of  target  motion  only  for  noise  limited  cases 
However,  to  detect  targets  in  the  presence  of  clutter,  target  motion  is  necessary  foi  the 
peak  to  migrate  into  the  clutter  free  region  of  the  profile  For  a  fixed  target  velocity,  a 
lower  PRF  results  in  further  migration  into  the  clutter  free  zone 

In  this  thesis,  a  novel  method  is  proposed  to  detect  moving  targets  in  clutter.  Targets 
shift  out  of  clutter  regions  due  to  range-doppler  coupling  The  resulting  shift  in  range 
confirms  the  presence  of  a  target  since  the  clutter  free  region  is  free  of  returns  from 
stationary  objects  Velocity  compensation  can  be  applied  :c  the  HRR  profile  to  determine 
the  true  range  of  an  extended  target  However,  there  are  two  problems  to  solve 
Specifically,  the  target  velocity  is  unknown,  and  the  portion  of  the  range  profile 
corresponding  to  the  true  range  is  obscured  by  clutter  These  problems  are  solved  as 
follows  The  cluttei  canceled  by  canceling  it  out  once  the  separation  of  clutter  and 
target  has  taken  place  Since  the  target  velocity  is  unknown,  compensation  can  be 
performed  for  several  different  velocities  sequentially  The  velocity  which  gives  the  Itighesi 
signal  magnitude  in  the  compensated  profile  is  the  correct  target  velocity 

Further  research  efforts  in  this  area  should  include  the  study  and  design  of  medium 
and  high  PRF  modes  of  operation  for  the  step  frequency  radar  In  addition,  the  properties 
of  the  step  frequency  waveform  using  the  ambiguity  function  need  to  be  studied  Finally, 
methods  to  further  refine  velocity  compensation  of  moving  targets  in  clutter  need  to  be 
developed 
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APPENDIX  A:  DESCRIPTION  OF  HRR  SIMULATION  PROGRAMS 


The  purpose  of  the  two  simulation  programs  HRRC  .M,  and  HRRN  .M  is  to  recreate  a 
synthetic  high  range  resolution  (HRR)  profile  in  clutter  limited  situations  and  noise  limited 
situations  respectively  HRRCAD.M,  the  program  which  simulates  clutter  cancellation 
followed  by  velocity  compensation  for  moving  targets  in  clutter  will  be  discussed  later 
Each  simulation  is  accessed  by  the  user  as  a  MATLAB  function  which  simulates  up  to 
three  targets  with  varying  ranges,  velocities  and  cross-sections  In  addition  to  the  range 
profile,  the  function  returns  the  processed  range  resolution,  and  the  unambiguous  range 
The  range  profile  is  displayed  as  a  plot  of  the  magnitude  response  of  the  inverse  FFT  in 
decibels  versus  the  bin  number  The  following  parameters  are  input  by  the  user  in  the 


order  shown  for  both  programs: 

In  “the  number  of  coherently  integrated  pulses 

2  fstep  =  the  frequency  increment  between  successive  pulses  (Hz) 

3  prf  =  the  pulse  repetition  frequency  of  the  radar  (Hz) 

4  rl  =  the  initial  range  to  the  first  target  (m) 

5  r2  =  the  initial  range  to  the  second  target  (m) 

6  r3  =  the  initial  range  to  the  third  target  (m) 

7  vl  =  the  radial  velocity  of  the  first  target  (♦  =  opening)  (-  =  closing)  (m/s) 

8  v2  =  the  radial  velocity  of  the  second  target  ’’  "  (m/s) 

9  v3  =  the  radial  velocity  of  the  third  target  "  "  (m/s) 


In  addition  to  these  direct  inputs,  the  following  parameters  are  easily  changed  within  the 
program  They  are  as  follows 
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I.  siir 


=  the  signal  to  noise  ratio  for  a  1  nr  target  (dB) 

2  scr  -  the  signal  to  clutter  j .  tio  fore  1  m*  target  (HRRC.M  only)  (dB> 

3  CO  =  the  compensaied  velocity  (+  =  opening)  (-  =  closing)  (,Ti/s) 

4  fo  =  the  carrier  frequency  (Hz) 

The  following  quantities  are  computed  and  presented  as  outputs: 

1  ru  =  the  unambiguous  range  (m) 

2  rstep  “  the  procesfcd  range  bin  size  (windowed)  (iii) 


The  output  parameters  can  be  changed  to  those  mere  relevant  to  the  user  such  as  the 
dispersion  factor  P,  and  the  range  bin  shift  factor  L 

Once  all  parameters  ate  specified,  the  first  step  of  the  program  creates  the  in  phase 
and  quadrature  samples  cf  the  target  return  signal  for  the  N  pulses  wnthin  the  burst  i  hese 
samples  are  written  in  complex  exponential  form  for  each  target  as 

s,(k)  =  A,fcxp(-j<p,(k))  (Al) 

where  A,  is  the  amplitude  of  the  return  from  the  kth  pulse  of  the  ith  target,  and  <j;, 
represents  the  corresponding  phase.  The  amplitude  of  each  target  return  is  assumed  to 
remain  constant  from  scan  to  scan  For  a  given  target,  the  phase  is  expiessed  as  follows 

(|)(k)  =  27tf(k)t(k)  (A2^ 

where  t(k)  is  the  carrier  frequency  of  the  kth  pulse,  and  t(k)  is  the  round  trip  time  for  the 
kth  target  return  The  carrier  frequency  of  the  kth  pulse  within  the  burst  is  given  by 

f(k)  =  f,.+(k-l)AF  (A  .3) 
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where  AF  represents  the  frequency  increment  between  successive  pulses  The  round  trip 
time  of  the  kth  target  return  is  given  by 


t(k)  = 


2.R(k) 

c 


(A.4) 


where  R(k)  is  the  target  range  for  the  kth  return  pulse  The  instantaneous  target  range  is 
specified  as 


R{k)  =  R,+(k-l)v(PRJ)  (A  5) 

where  v  is  the  radial  velocity  of  the  target  Substituting  the  expressions  for  the  target 
range,  the  target  round  trip  time,  and  the  carrier  frequency  into  equation  (A  2)  yields  tlie 
overall  expression  for  the  phase  of  each  return  from  a  given  target  Thus  the  phase  ())( k)  is 
rewritten  as 


(Kk)  -  — (f„  +(k  -  1)AF)(R„  +  (k  -  l)v(PRI))  (A  6) 

c 

A  velocity  compensation  factor  is  applied  to  each  target  return  signal  to  mitigate  the 
effects  of  the  range-doppler  coupling  caused  by  the  target's  radial  velocity  The  accuracy 
of  the  compensation  factor  depends  on  specifying  the  correct  target  radial  velocity  An 
uncompensated  profile  is  generated  by  setting  the  compensation  velocity  to  zero  The 
simulation  uses  the  same  compensation  factor  for  all  three  targets  Thus,  if  the  three 
targets  have  different  velocities,  none  of  the  targets  are  compensated  exactly  In  this  case. 
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the  best  possible  results  are  obtained  by  using  the  average  velocity  in  the  compensation 
factor  The  velocity  compensation  factor  is  given  as 


c 


(A7) 


where  co  indicates  the  compensation  velocity. 

Thus,  with  velocity  compensation,  the  kth  return  signal  for  a  given  target  becomes 

s,(k)  =  A,exp(-j4),(k))vc(k)  (A  8) 

In  both  simula'i'oiis,  the  SNR  represents  a  ratio  of  the  signal  power  to  the  noise  power 
as  shown  below  in  equation  (A  9) 


SNR  =  ^  (A  9) 

2a- 

The  reference  value  of  the  SNR  is  set  at  0  dB  for  the  1  m'  target  (T3)  prior  to  integration 
and  the  noise  source  is  modeled  as  white  gaussian  noise  with  a  variance  of  0  0002  Once 
the  noise  variance  is  set  and  the  SNR  is  specified  for  the  1  m*  target,  the  ainplitude 
corresponding  to  the  1  m‘  target  is  determined  from  equation  (A  9)  The  cross  sections  of 
the  remaining  two  targets  are  specified  as  10  m'  (T2),  and  100  m'  (Tl)  which  correspond 
to  target  SNR's  of  10  dB  and  20  dB  respectively  Since  the  target  cross  sections  are 
proportional  to  the  signal  power,  increasin  the  target  cross  section  by  a  factor  of  10  is 
equivalent  to  increasing  the  amplitude  A  Ly  a  factor  of  VlO  or  3  16 
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In  addition  to  the  actual  targets,  the  cluticr  simulation  (HRRC.M)  models  clutter 
by  including  an  additional  ten  closely  spaced  scatterers  within  the  anginal  range  bin  of 
width  ct/2  The  range  for  each  clutter  source  is  fixed  and  the  expression  for  a  single 
clutter  source  is  given  as  follows 


c,(k)  =  Cexp( 


-j47t(C-.(k-l)AF)R. 


) 


(A  10) 


where  C  represents  the  identical  magmtude  of  each  clutter  source.  The  value  of  C  is  based 
on  the  reference  value  of  the  signal  to  clutter  ratio  (SCR)  The  signal  to  clutter  ratio, 
(SCR)  represents  the  ratio  of  the  signal  power  of  a  1  m*  target  (T3)  and  a  single  clutter 
source  The  expression  for  the  SCR  is  given  as 


SCR- 


A*  2 
C'  2 


A* 

C' 


(All) 


where  C  is  the  identical  clutter  magnitude  for  each  clutter  source 

The  reference  value  of  SCR  is  set  at  -20  dB  for  the  1  m’  target  As  a  result,  the  signal 
to  clutter  ratio  for  the  10  m'  target  (T2)  and  the  100  m'  target  (Tl)  are  set  at  -10  dB  and 
0  dB  respectively  The  clutter  magnitude  C  is  determined  based  on  the  amplitude  A, 
computed  for  the  1  m’  target  (T3)  in  equation  (A  9),  and  the  specified  SCR  value  of  -20 
dB  for  that  target  The  clutter  magnitude  is  given  as 


C  = 


I  (aJ- 

\(SCR) 


(A  12) 
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The  overall  return  signal  for  the  kth  pulse  s(k)  which  consists  o>'  the  three  targets,  the 


noise,  and  the  clutter  (if  applicable)  within  one  original  range  bin  is  given  as: 


s(k)-^s,(k)  +  ]£c,(k)  +  an(k) 


(A  15) 


where  the  index  i  counts  the  number  of  targets  and  clutter  sources,  k  is  the  pulse  index, 
and  n(k)  represents  the  noise  source  whose  standard  deviation  o  equals  Vo'oobi  or  0.04 
Once  the  set  of  return  signals  are  generated  for  one  coherent  processing  interval,  the 
frequency  domain  data  are  then  multiplied  by  a  weighting  fiinction  to  reduce  the  sidelobe 
levels  of  the  inverse  FFT  A  Hamming  window  is  used  which  reduces  the  sidelobes  to  -43 
dB  below  the  main  lobe,  but  also  widens  the  main  lobe  by  2  dB  or  a  factor  of  1  33  As  a 
result,  the  processed  range  resolution  decreases  by  the  same  factor 

The  main  step  in  the  program  computes  the  Inverse  Fourier  Transform  of  the 
frequency  domain  samples  The  magnitude  of  the  inverse  FFT  represents  the  HRR  profile 
The  HRR  profile  is  plotted  with  the  IFFT  bin  number  on  the  x-axis,  and  the  magnitude  of 
the  IFFT  in  decibels  on  the  y-axis  The  final  step  computes  the  output  quantities  such  as 
the  unambiguous  range  and  the  processed  range  resolution  In  addition,  optional  output 
quantities  such  as  the  dispersion  factor,  and  the  bin  shift  factor  for  each  target  are  also 
determined  The  unambiguous  range  is  given  by  the  following  relation 


R.  = 


2AF 


(A  14) 
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where  AF  represents  the  freqy-ncy  increment  between  successive  pulses.  The  processed 
range  resoluuon,  is  given  as: 


Ar  = 


Kc 

2nAF 


(A  15) 


where  K  is  the  factor  by  which  the  mainlobe  is  widened  due  to  windowing  (K  =  133  for  a 
Hamming  window) 

In  contrast  to  the  two  programs  discussed  above,  HRRCAD  M  simulates  one  moving 
target  in  clutter  witliin  an  original  range  bin  of  width  ct/2  and  displays  the  output  in  tt  rms 
of  magnitude  rather  than  decibels  The  magnitude  pioi  is  chosen  to  better  illustrate  the 
resulting  improvement  in  the  range  resolution  The  input  parameters  to  this  program  are 


entered  in  the  following  order: 

1  n  =  the  number  of  coherently  integrated  pulses 

2  fstep  -  the  frequency  increment  between  successive  pulses  (Hz) 

3  prf  =  the  pulse  repetition  frequency  of  the  radar  (Hz) 

4  scr  =  the  signal  to  clutter  ratio  (d3) 

5  r  =  the  initial  range  to  the  target  (m) 

6  V  =  the  radial  velocity  of  the  target  (+  =  opening)  (-  -  closing)  (m/s) 


Once  the  initial  parameters  are  entered,  the  program  simulates  the  original  clutter 
limited  HRR  profile  using  the  same  general  methods  as  in  HRRC  M  The  simulation  then 
prompts  the  user  to  enter  the  highest  IFFT  bin  number  which  contains  clutter,  and  the 
compensation  velocity  for  the  target  These  two  parameters  are  referred  to  as  (bu)  and 
(col)  inside  the  program 
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The  high  bin  number,  for  the  clutter  region  determines  the  width  of  the  rectangular 
gating  function  whicfi  cancels  the  clutter.  The  gating  function  multiplies  IFFT  bin  numbers 
less  than  or  equal  to  the  high  bin  number  uith  zeros,  and  the  remainder  of  the  original 
HRR  profile  with  ones  to  implement  the  clutter  cancellation.  To  smooth  the  resulting 
sidelobes,  the  original  HRR  profile  is  then  multiplied  by  a  Hamming  window  prior  to  being 
transformed  once  again  to  frequency  domain  data  by  the  FFT 

Velocity  compensation  based  on  the  user's  estimate  is  applied  to  the  resulting 
fi  equency  domain  samples  The  compensated  data  are  then  transformed  into  the  final  HRR 
profile  by  the  inverse  FFT  The  resulting  profile  is  free  of  clutter,  and  indicates  an  increase 
in  the  target's  signal  strength  provided  that  the  correct  compensation  velocity  has  been 
entered  In  addition,  the  program  returns  the  unambiguous  range  R,,  and  the  processed 
range  resolution  Ar 
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APPENDIX  B:  HRR  SIMULATION  PROGRAM  CODES 


function  (ru.rstepl  =  hrm(n,fst9p,prf,r1,r2.r3.v1,v2.v3) 

%  This  program  is  a  MATLAB  function  which  computes  the  HRR  profile  of 
%  a  single  range  bin  in  noise  limited  situations 
%  The  following  data  are  inputs  for  each  of  tht  3  targets: 

%  Number  of  pulses  (n) 

%  Frequency  step  size  (fstep) 

%  Pulse  Repetition  Frequency  (prf) 

%  Initiai  Target  Ranges  (r1 .  r2,  r3] 

%  'arget  Radial  Velocity  (vl .  v2,  v3)(-»  opening)(-  closing) 

%  The  following  parameters  can  be  changed  within  the  function 
%  Signal  to  Noise  Ratio  in  (dB)  for  each  target  (snr) 

%  Nominal  carrier  frequency  (fo) 

°/b  Velocity  matched  for  compensation  (co) 

%  The  outputs  are  Ru  =  c/2*fstep,  and  Rstep  =  Kc/2*n*fstep, 

%  which  are  the  unambiguous  range  and  processed  range  bin  size. 

%  The  FFT  size  is  equal  to  n 

%  The  frequency  domain  data  are  weighted  by  a  Hamming  window.  (K=1.33) 

%  The  noise  is  white  gaussian  noise  with  zero  mean  and  variance  equal  to  0002. 

%  Initialize  parameters 

co  =  0; 

fo  =  Ie9: 

snr  =  0; 

x=[0;ri-11, 

c  =  3e8, 

K  =  1.33; 

PRI  =  1/prf, 

%  Convert  the  SNR  from  decibels  to  magnitude 
a  =  10''(snr/10): 

A3  =  sqrt(2*a*.0002); 

A2  =  sqrt(in)*A3, 

A1  =  sqrt(1  0)*A3: 

%  Generate  the  frequency  domain  data 
for  k=1  n, 

f(k)  =  fo  +  (k-1)*tstep; 

vc(k)  =  expO*4*pi*co*f(k)*(k-1)*PRI/c); 

Sl(k)  =  Arexp(-(4*pi/c)*f(k)*(r1  *(k-1)*vrPRI))*vc(k), 
s2(k)  =  A2*exp(-(4*pi/c)*f(k)*(r2  ♦(k-1)*v2*PRI))*vc(k) 
s3(k)  =  A3‘exp(-(4*pi/c)*f(k)*(r3  •►(k-1)*v3*PRl))‘vc(k); 
s(k)  =  (sl(k)+s2(k)+s3(k);  +  sqrl(.0002)*randn; 
end 
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%  Compute  the  IFIFT 
h  =  hamming(n); 
sw  =  s  *  h'; 

S  =  ifft(sw); 

mag  =  0bs(S); 

dbmag  =  20*log10(mag); 

%  Compute  output  quantities 

ru  =  c/(2*fstep); 

rstep  =  (K*c)/(2*n*fsiep): 

P  =  v1*n*PR|/rstep; 

B  =  n*fstep; 
fc  =  fo  +  B/2; 

L  =  (fc/B)*P; 

®/b  Plot  the  HRR  Profile 

plot(x,dbmag) 

grid 

xIabelCIFFT  bin  number  ') 

yiabel('Uagniiude  response  in  dB') 

gtextCTi*) 

gtextCT2') 

gtextCTS') 


function  (ru.rstep]  =  htTc(n,fstep,prf,r1.r2,r3,v1.v2.v3) 

%  This  program  is  a  MATLAB  function  which 
%  computes  the  HRR  profile  in  clutter  limited  situations 
%  The  following  data  are  inputs  for  each  of  the  3  targets; 

%  Number  of  pulses  (n) 

%  Frequei;cy  step  size  (fstep) 

%  Pulse  Repetition  Frequency  (pd) 

%  Initial  Ranges  (r1 .  r2.  r3) 

%  Target  Radial  Velocity  (v1,  v2,  v3)(+  opening)(-  closing) 

%  The  following  parameters  can  be  changed  within  the  function 
%  Nominal  earner  frequency  (fo) 

%  Signal  to  Noise  Ratio  in  (dB)  for  each  target  (snr) 

%  Signal  Clutter  Ratio  in  (dB)  for  each  target  (scr) 

%  Compensated  velocity  (co) 

%  The  outputs  are  Ru  -  &'2*fstep.  and  Rstep  =  Kc/2*n*fstep. 

%  which  are  the  unambiguous  range  and  processed  range  bit.  size 
%  The  frequency  domain  data  are  weighted  by  a  Hamming  windov/.  (K  =  1  33) 
%  The  noise  is  ^ite  gaussian  noise  with  zero  mean  and  unit  variance 

%  Initialize  parameters 
fo  =  1e9. 
snr  =  0; 
scr  =  -20: 

CO  =  0: 

K=  1.33: 
c  =  3eS: 
x=  (0:n-11; 
a  -  10  ^(snr/10); 

A3  -  sqrt(2‘a*.0002): 

A2  =  sqrt(10)*A3; 

A1  =  sqn(100)*A3; 
b=  10 ''(scr/10); 

C  sqrt((A3V)/(b)); 

%  Generate  the  frequency  domain  data 
for  k=1:n, 

f(k)  =fo-*- (k-1)*fstep; 

vc(k)  =  exp(j*4*pi*co*f(k)*(k-1)*PRI/c): 

«;1(k)  =  Arexp(-(j*4*pi/c)*f(k)*(r1  +(k-1)-v1*PRI))*vc(k); 

s2(k)  =  A2*exp(-a*4*pi/c)*f(k)’(r2  +(k.1)*v2*PRI))*vC(k); 

s3(k)  =  A3'‘exp(-(j*4*pi/c)*f(k)*(r3  ••■(k-1)*v3*PRi))*vc(k). 

c1(k)  =  C'exp(-0"4*pi/c)*f(k)*1), 

c2(k)  C*exp(-0*4‘pi/c)*f(k)*2). 

c3(k)  =  C*exp(-0*4*pi/c)*f(k)*3): 

c4(k)  =  C*exp(-(j*^‘pi/c)*f(k)’5), 

c5(k)  C*exp(-(i*4*pi/c)*f(k)*6): 

c6(k)  =  C*exp(-(j*4'*pi/c)*f(k)*7): 

c7(k)  =  C*exp(-(j*4‘pi/c)*f(k)*9); 

c8(k)  -  C*exp(-0*4*pi/c)*f(k)”11) 

c9(k)  =  C*exp(-0*4*pi/c)*f(k)*13): 

c10(k)  =  C*  'xp(-G*4*pi/c)*f(k)*15); 
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s(k)  =  (s1(k)+s2(!i)+s3(k))-*-(c1(k)+c2(k)+c3(k)+c4(k)-*-c5(k)+c6(k)  ... 

+  c7(k)>c8(k)+c9(k)+c10(k))  sqrt(  0CO2)*randn: 


end 

%  Compute  the  IFFT 
h  =  hamming(n); 
sw  =  s  .*  h'; 
mag  =  abs(S): 
dbir.ag  =  20*log10(mag). 

%  Compute  the  output  quantities 

ru  -  c/(2*fstep); 

rstep  =  K*c/(2*n*fstep); 

pi  =  v1*n*PRI/r5tep; 

p2  =  v2*n*PRI/rstep: 

p3  =  v3*n*PRI/rstep: 

P  =  [p1  p2  p3); 

B  =  n*fstep; 
fc  =  fo  *  B/2: 

L  =  (fc/B)  .*  P; 


%  Plot  the  HRR  Profile 

plot(x,dbmag) 

grid 

xlabelC'FFT  bin  number  ') 
ylabelCMagnitude  response  in  dB') 
gtextv'TI) 
giexi(T2') 

9lext(T3') 
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function  (ru.rstep]  =  hrTcad(n.fsiep.prf,scr,r.v) 

%  This  program  is  a  MATLAB  function  which  computes 
%  the  HRR  profile  of  a  single  moving  target  in  clutter 
%  This  program  includes  clutter  cancoilation 
%  followed  by  velocity  compensation. 

%  The  input  data  consists  of  the  following: 

%  Number  of  pulses  (n) 

%  Frequency  step  size  (fstep) 

%  Pulse  Repetition  Frequency  (prf) 

%  Signal  to  Clutter  Ratio  for  1  square  meter  target  (scr) 

%  Initial  Target  Range  (r) 

%  Target  Radial  Velocity  (v)  (+  opening)  (-  closing) 

%  The  following  parameters  can  be  changed  within  the  function 
%  Amplitude  of  return  signal  for  1  square  meter  target  (A1) 

%  Nominal  carrier  frequency  (fo) 

%  The  following  parameters  are  input  by  the  user  after  the 
%  original  HRR  profile  is  displayed: 

%  Highest  IFFT  bin  containing  clutter  (bu) 

%  Compensation  velocity  (col) 

%  The  ou*iuts  are  Ru  =  c/2*fstep,  and  Rstep  =  K*c/2*n*fslep. 

%  which  a  a  the  unambiguous  range  and  processed  range  bin  size. 
%  The  FFT  size  is  equal  to  n 

%  initialize  parameters 
fo  =  1e9i 

V*  *«  1 ' 

^  |W...  .j, 

c  =  3e8, 

PRI  =  1/p.-f; 

A1  =  1: 

b=  10  "'(scr/IO), 

%  B  is  the  clutter  magnitude 
8-  sqrt((A1''2)/b) 

for  k=1:n, 

f(k)  =  fo  +  (k-1)*fstep, 

s1(k)  =  Arexp{-(j*4-pi/c)*f(k;*(r1  +(k-1)*vrPRI)); 

c1(k)  =  B*exp(-(j'4'pi/c)*f(k)*1): 

c2(k)  =  B*exp(-(j*4*pi/c)*f(k)*2); 

c3(k)  =  B*exp(-(j*4*pi/c)*f(k)*3): 

c4(k)  =  B*exp(-0*4*pi/c)'f(k)*5), 

c5(k)  =  B*exp(-(.i*4*pi/c)*f(k)''6): 

c6(k)  =  B*exp(-0*4*pi/c)*f(k)'*7), 

c7(k)  =  B*exp(-(j*4*pi/c)*f(k)*9), 

c8(k)  =  B*exp(-(j*4*pi/c)*f(k)*11). 

c9(k)  =  B'exp(-rj*4*p:/c)*f(k)*13). 

c10(k)  =  B*exp(-G*4*pi/c)*f(k)*15), 

£(k)  =  s1(k)  +  (c1(k)+c2(k)+c3(k)+c4(k)+c5(k)'*-c6(k) 
+c7(k)+c8(k)'«-c9(k)+c1 0(k)): 
end 
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%  Compute  original  range  pronie 

%  Determine  output  quantities,  Ru  and  rstep 

S  =  iffl(s); 

mag  =  abs(S); 

ru  =  c/(2*fstep); 

rstep  *  K*c/(2*nfstep), 

plot(x.mag) 

grid 

titieCOriginal  Range  Profile’) 
xlabelCIFFT  bin  number  (origins!) ') 
ylabelCMagnitude  response  *) 
gtextCTt') 

%  Obtain  user  input  for  clutter  cancellation  and 
%  compensation  velocity 
bu  =  inputCEnter  high  bin  number  for  clutter  ’) 
col  =  inputCEnter  velocity  for  target  1') 
ssl  =  |zeros(1:bu)  ones(1:n'bu}]  .*  S; 

%  Wi  ^dow  the  range  domain  data 
h  =  hamming(n): 
ssl  =  ssl  .*  h': 

for  k=1  n, 

f(k)  =  fo  +  (k-l)*fstep; 
vc1(k)  =  expG*4*pi*co1*f(k)*(k-1)*PRI/c); 
er>rj 

%  Apply  velocity  compensation  to  transfomied  data 
SS1  =fft(ss1)  *vc.1; 

%  Compute  final  range  profile 
FSS  =  itfKSSI); 
finpiot  5=  (abs(FSS)); 
figure 

plot(x,finplot); 

titleCCompensated  HRR  profile') 
xlabeK'IFFT  Bin  No’) 
ylabelCMagnitude  response’) 
grid 

gtext(’Tl’) 
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